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ABSTRACT 


Interfaces between particles in particulate materials produce a microstrucrure 
consisting of line segments of varying lengths on the plane of polish. The interfaces 
between particles can have variuous sizes and shapes. In this work a stereological analysis 
has been conducted on circular and elliptical interfaces. Relationships between some 
stereological parameters, such as, mean tangent diameter and mean size of interfaces, and 
the length of segments have been derived. Stereological methods have been developed to 
unfold the shape and size distribution of inter-particle interfaces from the distribution of 
line lengths ( which can be measured from observed microstructures ). 

In order to test the validity of the above analysis microstructures were generated 
by simulating discs of circular or elliptical shapes (discs represent the interfaces between 
particles) and various sizes distributed in 3-D space. This 3-D space was then sectioned 
by a random plane (which simulates the plane of polish) to obtain a microstructure 
consisting of line segments (intersections between discs and the random plane). Using the 
above methods, the line length distribution obtained from these microstructures was 
analyzed to extract information on the shape and size distribution of discs. 

From the above it has been concluded that the stereological relationships and the 
methods developed for the transformation of line length distribution to disc size 
distribution are valid. It has also been shown that it is also possible to shed some light on 
the shape of the interfaces. 


CHAPTER 1 

INTRODUCTION 


The structure of materials consists of several components such as grain-boun dar ies, interfaces, 
precipitates etc. The geometrical characteristics, such as size, shape etc., of these components have 
profound influence on mechanical properties. One of the well established correlation between 
microstructure and properties is the Hall-Petch [1] relationship between strength and grain size of 
metals and alloys. Size and size distribution of particles and grains have also been shown by many 
researchers to have an influence on various mechanical properties. 

Figure 1.1 shows an example of a microstructure in a fully annealed polycrystalline specimen of 
3 % silicon iron. The sample shows many grain-boundaries or interfaces between the grains. These 
interfaces are planes in 3-dimensional structure and have different shapes and sizes. Figure 1.2 
illustrates an SEM micrograph of a Ti-Al alloy. In the microstructure the particles of TiAl 3 have 
precipitated out from the alloy. The 2-D profile of the particles visible in the plane of polish 1 ' , shown in 
this figure are actually sections of 3-D polyhedra. 



* A plane which intersects a sample to reveal its 2-D information, is called a plane of polish. 



..V 



Figure 1 .2 . Fully annealed polvcrystalline specimen of 3% silicon iron, etched to reveal 
grain-boundaries. 


In powder processing, powder is compacted and sintered for manufacturing purpose. At the 
beginning of compaction, the powder has a density approximately equal to apparent density. For a 
loose powder there is an excess of void space, no strength and a low coordination number. As the 
pressure is applied, the first response is rearrangement of the particles with filling of large pores, giving 
a higher packing coordination. Increasing pressure provides better packing and leads to decreasing 
porosity with the formation of new particle interfaces contacts [ 3 ]. During the initial stage of sintering 
a rapid growth of the interparticle necks or bonds is observed. In the intermediate stage, the pore 
structure become smoother and has interconnected cylindrical necks. In the latter stage of sintering, 
grain growth occurs, giving a larger average grain size with fewer grains. Figure 1.3 shows the 
progressive densification and grain growth during sintering. The interfaces between powder particles 
formed in the initial stage of sintering are converted into grain-boundaries in the final stage of sintering. 
In the following, the formation of snow as an example of a particulate material is briefly discussed. 

The process of compaction, sintering and subsequent densification occurs in nature during the 
formation of snow. Snow is basically a granular material consisting of ice particles which are bonded 
to their immediate neighbours through inter-particle bonds. Pile-up of snow on the ground undergoes 
sintering under pressure generally at temperatures close to its melting point. The ceformaxion 
behaviour of snow on application of pressure can be described in the following manner: 




c) final 


Figure 1.3 : The progressive densification and grain growth during sintering. 


Bonds established between coordinating particles in a loose pack of snow are not capable of 
withstanding high loads. Many bonds between ice particles, therefore, undergo an immediate fracture 
on the application of external load. Fracture of a majority of initial bonds provide an opportunity for 
more efficient packing of panicles by their readjustment and the mass now occupies a smaller volume. 
The extent of bond fracture and panicle readjustments, however, depends on the initial panicle 
arrangement and the magnitude of the externally applied load. Fracture of bonds and readjustment is 
called stage I of deformation of snow. 



From stage II onwards the structure deforms without any further readjustment of its particles. 
Under these conditions, the irreversible deformation of snow continues to occur bv the time dependent 
growth of inter-particle bonds. The strain rate in stage II keeps decreasing as the deformation proceeds 
because the effective pressure on inter-particle bonds decreases as their total area on a given particle 
increases. Soon the strain rate reaches a steady state and at this point, stage II of deformation 
terminates. Figure 1.4 shows the vertical snow profile of microstructure immediately after mechanical 
processing. The ice particles are shown in white while the pores appear black. The pore volume can be 
reduced by further applying pressure or by increasing the time of mechanical processing [ 5 ]. 



Figure 1 .4 : Verticle snow profile of microstructure after mechanical processing, 
(a) binary image of microstructure (b) grain skeleton 


In the formation of snow from ice particles, Kry [ 7 ] assumed the bonds between ice particles 
to have the shape of circular discs. In snow structure, it is necessary' to find out the bond radius 
between two grains because during deformation the largest stress concentrations occur at the bond 
since the cross-sectional area is smaller than that of associated grains. Owing to the increased stress 
concentration through the neck region, any plastic flow which occurs will generally occur in the necks 
rather than in the grain bodies. These necks, therefore, play a major role in determining the behavior 
during loading. The changes in the neck sizes can be used to measure the plastic flow. 


Another significant factor controlling the mechanical behavior of snow is the number of bonds 
per grain ( 3-dimensional coordination number ). As the number of bonds per grain increases, more 
load can be supported, thereby providing a mechanism for reducing the stress concentration in a given 
bond. If, for a given loading, the stress in a bond can be reduced, then the probability of its breaking 
will be reduced. If the stress is reduced in a large enough number of bonds, the overall load-bearing 
capacity of the snow can be increased [ 8 j.The degree of bonding between grains and the pore length 
(or the mean free distance between grains ) is considered to control the mechanical properties. 
Connected with these are the number of grains and bonds per unit volume as well as the free grain 
surface area per unit volume. Larger the number of bonds and/or bigger the radius, the higher the 
resistance the material offers to the deformation and higher will thus be the modulus as shown in figure 
1.5. 



Figure 1.5: Dependence of Young’s modulus on free grain-surface area per unit 
volume. 


The mean 3-dimensional bond radius is a significant parameter for pressure sintering as for 
bond strength and fracture. The bond radius comes out to be significant parameter for the deformation 

to 

of snow. However, It turns out that instead of individual bond areas formed on a single particle, it is the 


ratio of total bond area and the initial free surface area of a grain which plays a critical role in affecting 
the deformation of snow. The variation of this ratio i.e. St/S 0 , during deformation of snow has been 
shown in figure 1 . 6 .. 
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Figure 1.6: Sf. S 0 versus time at varying pressures. 

This study attempts to analyze the size and shape of interfaces between particles, grains etc. 
and also thin discs like precipitates in multi-phase and particulate materials. The present study focuses 
on the shape and size analysis of interfaces/grain-boundaries in particulate materials but is restricted 
only to circular and elliptical shapes. In this study, the interfaces between particles appear as lines when 
such a 3-D structure is intersected by a sectioning plane (i.e., plane of polish), these lines have varying 
lengths as can be observed in figure 3.1. From the distribution of these line lengths it is possible to find 
the shape and distribution of interfaces/'grain-boundaries in 3-dimensional structure by taking various 
assumptions into account. 

In this work, an attempt is made to develop analyses techniques for the line length distribution 
to obtain the following objectives. 

1 . To develop relationships for some basic stereological parameters, such as, surface area / volume 
(S v ), Number of interfaces / volume (NY). Number of intercepts / length (NO, Number of profiles 
/ area (N A ), Length of profile / area (L A ) and average value of fee lengths. 
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2. To develop a transform for estimating the distribution of circular and elliptical disc sizes in the 3-D 
structure from the measurement of the line lengths distribution. 

3. To analyze the shape of the discs or the interfaces. 

In the next chapter, some basic principles of stereology which have been used in subsequent 
analyses, are discussed. 


CHAPTER 2 


FUNDAMENTAL CONCEPTS IN STEREOLOGY 


Stereology is a body of mathematical methods relating 3-dimensional parameters defining the 
structure to 2-dimensional measurements obtainable on sections of the structure. Using stereology 
it is possible to give quantitative interpretation of 2-dimensional obse nations in terms of 3- 
dimensional observations. Mathematical methods have been developed which allow obsenations 
made on sections to be interpreted in terms of spatial structure of the specimen [11]. 

This chapter deals with some of the fundamental concepts in stereology. which are discussed 
in the subsequent sections. 


2.1 Random Sectioning 


Figure 2. 1 shows a specimen intersected by a plane of polish. The plane of polish is a plane 
which cuts the specimen to provide a surface for microstructural obsen'ations. The orientation 
(0,cj)), of the sectioning plane is defined by its normal which makes an angle 0 w ith z axis and whose 
projection on XY plane makes an angle 0 with x-axis. 

Consider a line OA which is normal to the section plane and passes through the origin as 
shown in Figure 2.1. The section plane T is defined by the distance from the origin ai which the 
plane makes an intersection with the specimen. The distance from the origin to the section plane in 
the direction (0,0) is given by r with range {-«< r < *>). Thus the plane T is uniquely defined by T 
= T(0,0.r). Here 0 and 0 have the range {0<9<~. 0<o< 2it}. In order to find a random plane 





which intersects the specimen S, assign random values to 0,0 and r within their specified ranses. By 
assigning random values to 0,0 and r, one can obtain 'Isotropic Uniform Random' (IUR) planes 
[ 12 ]. 


2.2 Averaging with Respect to Orientation 


Consider a sphere of unit radius with centre at the origin of XYZ coordinate, as shown in 
figure 2.2. This figure shows three rays 1, 2 and 3 which are emanating from the origin. For random 
orientation the intersecting points between the rays and the sphere must be homogeneously 
dispersed over the sphere surface. The sphere surface can be subdivided into a large number of very 
small units of area 'ds'. Now, the probability that a ray crosses one of these units 'ds' can be written 
as 



ds 

Surface, area, of . sphere 


ds 

4tc 


( 2 . 1 ) 


Figure 2.3 shows a line whose orientation is defined by two angles Band 0, where 9 is the . 
angle between the line and z-axis and b is the angle between the x-axis and the projection of the line 
onto x-y plane. The range of variation of angles 9 and 0 is 0 to it and 0 to 2k respectively. In order 
to find out the probability that a line will have orientation (0,0), i.e. it will lie in the range { 0/(9+d0); 
0/(0+dp)}, a small quadrangular on the sphere surface is made by sweeping angles a0 and dp 
around the line OA in figure 2.3. The area of this quadrangular at a particular orientation (0,0) is 
given by 

(sinQ.dQ ).d§ (2.2) 


The area of this small quadrangular is not equal but is different for different orientations. This 
wedge-shaped region tapers to width 0 as the line approaches to z-axis. Le. as 0— >0. 




Probability of a line oriented in the range {0/(B-fd8); <i)/((j)4d(j)) } can now be defined as the 
ratio of the area of the small quadrangular shape region to the area of the sphere and is given bv 


Pr(Q,<$) 


sinQ.dQ.dty 

4k 


(2.3) 


The random orientations are uniformly distributed over the orientation sphere with probability 
density given by equation (2.1). 


In order to find the average value of any parameter pm(9,4>) which varies with orientation, 
integrate the product of parameter and probability over the whole range of 0 and p. Thus, the average 
value of parameter Pm (0,<j>) is given by 

P m = (2.4) 

4k 


Averaging with respect to orientation in 2-dimension follows a procedure which is analogous 
to that in 3-dimension. Consider a circle of unit radius with centre at the origin of XY coordinates 
as shown in figure 2.4. The orientation of a line in 2-dimension is defined by only one angle 9 which 
the line makes with x-axis. as shown in figure 2.4. Probability of finding a line in the range 
{9/(9+d9)}is given by 

dQ 

P(9) = (2.5) 

2k 

The mean value of any parameter p m (0) is given by the integration of the product of 
parameter and probability P(9) over the whole range of 0, i.e., 0 to 2 tc : 

p = —C-’p JQ)dB 
Pm 2lt ’ 


(2.6) 



2.3 General Conditions for the Validity of Stereological Principles: 


In order to choose a random point from the sample S from which the random plane passes, it 
is necessary to identify a sufficiently large number of uniformly distributed points in S. Then, one of 
these points is selected randomly. It is also necessary to ensure that the section plane is randomly 
oriented i.e. it has random values of 9, p and r. 


Very often, however, the specimen S is itself a sample from a much larger body W and the 
area of the section used for microscopic study is very small compared to the size of W. It is called as 
"extended case". For the extended case, where the section is cut from an IUR random sample of 
the very large specimen, arbitrary sections yield unbiased ratio estimators. 


The analogous conditions prevail with respect to "cutting" the specimen with a test line. In 
the extended case arbitrary line sections of the sample S out of the specimen W yield unbiased ratio 
estimators, provided S is an IUR sample of W [ 14 ]. 


2.4 Basic Stereological Parameters 


For every structure, there are two types of parameters. The parameters such as surface area 
of second phase particles per unit volume and number of particles per unit volume, which describe 
the properties of 3-D structures . While other parameters such as number of intercepts per unit 
length and length per unit area' termed as 2-D parameters, describe the properties of microstmcture 
observed on the plane of polish. There is a certain relationship between the 2-D and 3-D 
parameters. Some of the basic relations are discussed in the following subsections. 


2.4.1 Surface density, Sv 



Consider a cube of dimension 'I. The cube contains second phase particles, as shown in 
Figure 2.5. The second phase particles are randomly oriented in the cube. Suppose the surface of 
the second phase particles is represented by a polyhedral surface formed of very small polygons of 
equal area ‘ds’, as shown in figure 2.6. If ‘N : is the number of polygons covering the surface then 
the total surface area S T of the particle is given by 

St = A Ids (2.7) 


The surface of the second phase particle is assumed to be isotropic, i.e. that the surface elements 
‘ds’ takes all possible orientations in space. Choose randomly one of these surface elements and 
suspend it in the cube as shown in figure 2.7 . Section this structure with a test line parallel to z- 
axis. The probability that a random test line parallel to z-axis will intersect the surface element ‘ds’ 
will be given by the ratio of the area of the orthogonal projection of ‘ds’, denoted as Ap, to the 
orthogonal projection of cube, which is 1 : : 




A p (Q,§ )_ 
f ' 


( 2 . 8 ) 


and the average probability of intersection is given by 


PAD 



(2.9) 


where A is the projected area averaged over all the orientation. 


It is clear from figure 2.7 that the projected area Ap(0,b) is independent of angle o and 
depends only on the angle 0. therefore, the value of Ap(0,o) is given by 

A p (Q,§) = ds.cosQ 


( 2 . 10 ) 



i 


Figure 2.5 : Surface density estimation from intersection count on test lines 

[15]. 


Sfa) $(o)*N'dJ 



] Figure 2.6 : Polyhedral approximation of smooth surface [15]. 
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Figure 2.7 : Probability of hitting a surface element and its dependence on 
9 and 4> [15]. 
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The average value of A p (0) is given by the integration of the product of A p (0.o) with 
the orientation probability density P^b) over the entire range of 0 and 6 . Thus, from 
equations (2.4) and (2.10), the average projected area is given by 

; 1 2® it J 

A = — f . J ds.cos.Q.sinQ.dQ. do = — ds (2.11) 

' 471 0 0 '2 

From equations (2.9) and (2.11) the probability of intersection between the surface 
element and one of the test lines is given by 

„ 1 ds 

PAD = Y7 (2A2) 

In the case of n test lines each of length L the expected number of intersections E(I), 
between test lines and particle becomes 

E(I) = P (I). N.n = —~-N.ii =-S — 

2 l 2 2 r 

or, 

E(I) _ 1 S 

L ~ 2V 

From above 

S v =2.N l (2.13) 

where S v is the surface area of second phase particles per unit volume = S/V and N L is the 
number of intercepts made by test lines per unit length = E(I)/L . 

2.4.2 Boundary Density L A 

In the fisure 2.8a profile boundaries form a curve of length L, 
if the curve is smooth it can be cut into infinitesimal straight line segments of equal length 'dL 


such that 



(2.14) 


L = \..dL =n. dL 

Over. the. curve 

where, n is the number of infinitesimal line segments. The curve is assumed to be isotropic 
such that all the orientations of the segment dL are equally likely. A grid of parallel 
equidistant test lines are placed onto the section plane, as shown in figure 2.8b. The distance 
between the test lines (d) is kept greater than the length of line segment in order to avoid 
more than one intersections of a line segment with test lines. Suppose N is the number of 
intersections made by line segments on to the test lines. The probability of intersection of a 
line segment with a test line Ti is given by 

p r {dL/ Tij - N/n (2.15) 

The value of m, the projection of dL onto the orthogonal complement to the grid is given 
by 

m = dL cosQ (2.16) 

The probability that the element dL intersects a test line Ti of the grid is given by 

, , dL.cosQ 

p r fdL/ Til = m/d = (2.17) 

d 


For isotropic orientation of dL, 0 must vary uniformly over the range {0 < 8 < 2 ji} 
and the average probability of intersection is given by (using equation (2.6) ) 


PjdL/T,} = 2-J TPrldL/TJdB 

2K 


(2.18) 


From equations (2.15), (2.17) and (2.18): 


2.dL 

%.d 


N/n 


(2.19) 


From equations (2. 14) and (2. 19), the length. L of the boundary in figure 2.8a is given by 


L = -N.d 
2 


( 2 . 20 ) 




If both the curve and the test grid Ti are enclosed in an area A, the length of test line is : Li 
= A/d. Substituting the value of d in equation(2.20) the following relationship between the 
curve length density, La. and the number of intersections per unit test line length, Nl. is 
obtained as: 


La 


L _ n N 
A ~ ~2'T, 



( 2 . 21 ) 


2.4.3 Relationship between Surface Density in 3-D Space and Length 
Density in 2-D Space. 

It is clear that on combining equations (2.13) and (2.21), one gets the following 
relationships between length density L A and surface density Sv : 

Sv = ~La (2-22) 

K 

Stereological relations given by equations (2.13), (2.21) and (2.22) are statistically exact 
relationships. These relationships are orientation independent and consequently do not depend 
on the shape of the objects. However, there are other stereological parameters which can be 
estimated only if the shape of the objects are known. Some of these parameters are discussed in 
the following sections. 

2.5 Tangent Diameter 

"Tangent diameter of a particle in a particular direction is the distance between two 
parallel planes having their normals in the same direction and are tangent to the particle , as 
shown in figure 2.9. Suppose Q is a particle whose centre coincides with the origin of the 
XYZ coordinate. The distance between two parallel planes tangent to the particle is denoted 
as the tangent diameter H. as shown in the figure 2.9. The notation for tangent diameter m a 

particular direction (0,0) is given by H(0,<j>) . 



2.5.1 Average Tangent Diameter 


Average tangent diameter means the average of all the tangent diameters found along 

different directions (6.6) . Therefore, the average tangent diameter, H is given by (using 
equation (2.4) ): 

H = m 6,<t> J sinB dd (2.23) 

4tc 

It is clear from above discussion that the average tangent diameter is dependent on the 
shape of the particle. 

The value of average tangent diameter for some particles are given below [16] 

Table 2.1: Average Tangent Diameters, H 


Particle Shape 

Characteristic Dimension 

Average Tangent Diameter 

Sphere 

radius R 

2 R 

Cylinder 

length 1, radius r 

(1 + pR)/2 

Hemisphere 

radius R 

(1 +p/4)R 

Parallelepiped 

edge lengths: a.b,c 

(a + b + c)/2 

Cube 

edge length: a 

3 a/2 


2.6 Numerical Density of Particles 


In this section the relationship between the number of profiles seen on a random 
section and the number of objects can be established [ 17 ]. The objects in this case are 
discrete particles which are contained in the unit volume of the structure, as shown in the 
figure 2. 10 . Assume a large number of particles (no) to be enclosed in a very large space 
(V T ) such that the numerical density of the particles(N v ) is given by 

Nv = Ho/Vt (2.24) 

As shown in figure 2. 10, the sample is sectioned by a test plane of area At- On this n p 
number of particle profiles are observed which have their centers within At- The numerical 
density of profiles N A is given by 

N a = 7i,/ Ar (2-25) 

In order to find a relationship between N A and Nv, the following conditions need to be 
specified. 

1 . The particles must be convex solids, so that a plane intersecting the particle at 
any orientation will produce one profile only. 

2. The particles must be randomly oriented within the space and their number 
must be so large that all orientations between particle and test plane are equally likely. 

3. The size of the sample space should be such that the number of particles it 
contains is very large. 

The sample space is randomly sectioned by a plane parallel to X-Z plane and is 
defined by a cube of volume V T = l 3 , so that the section Area is given by 

= l 2 


At 


(2.26) 




Figure 2.10 : 


Sectioning of randomly oriented particles by a section plane 
[18]. 



Figure 2. 1 1 : Probability of cutting a random particle [1 8]. 



Probability of cutting a single particle Y by the random section plane is given by the ratio of 
the projections of both particles and containing space on to the orthogonal complement to the 
section plane as shown in figure 2. 1 1 . 


PAY /TJ = H(Y)/l (2.27) 

Considering all possible orientations the probability to cut any one of the particles is given bv 

P r {P article / Jj - H/l (2.28) 


The expected number of profiles on section Ti can be written as 

n p =n 0 H / l 


Dividing both sides by l 2 


n p / l 2 — n 0 H / Z 3 


(2.29) 


(2.30) 


Numerical density, N A . of particles on 2-D section can be related to the numerical density of 
particles, Nv, in 3-D space using equations (2.24), (2.25) and (2.30) : 

N a =N v .H (2.31) 


If the particle population represents a mixture of different sizes and shapes of particles, one 
can specify J{j for ea ch size and shape class and the relative frequency fj for each class. The total 
numerical density N v is the sum of the numerical densities of the particles belonging to each class 

(N vj ). 

N v =lN v (2.32) 

j 1 

For each class the numerical profile density is given by 

N Ai = N v . .~H i 


(2.33) 



And the total profile density can be written as 


N,=2N A =1N V .H, (2.34) 

j J j } 
or ? 

N A =N v .[Zf t .H l ] (2.35) 

j 

where, / ; . = N v / N v , is the frequency of the j* class of particles. The term in the bracket of 
equation (2.35) represents the weighted average of the mean tangent diameters. 


2.7 Size Distribution of Section Profiles when Spheres are Randomly 
Sectioned 


The intersection of a sphere with a plane produces a circle on the section plane. The size of 
the profile (circle) depends on the distance of plane from the centre of sphere. The maximum size of j 
profile is obtained when the plane passes through the centre of sphere. For a sphere of radius R s .the 
profile size can have radius between 0 and R, depending on the distance of the plane from the centre 

j 

of sphere. 

2.7.1 Spheres of Constant Size 

Figure 2. 12 shows a sphere of radius Rs placed in the XYZ coordinate system with its centre j 

at the origin. A plane perpendicular to the z-axis makes an intersection with this sphere. Each | 

section will produce a circular profile whose radius r will depend on the position of the section plane 

with respect to Z-axis . The value of radius r of a circular profile obtained by intersecting a sphere of 

...... ... ... | 

radius R s at a height z from the origin by a section plane is given by j 






Differentiating the above equation with respect to r gives: 



(2.36) 


(2.37) 


Probability of cutting a sphere of radius R s between height z to z+dz by a random horizontal 
section plane is given by 

P( Z ) = — (2.38) 

2R s 

where 2R S is the tangent diameter of the sphere of radius R s as shown in Figure 2.12. The radii of 
circular profiles, obtained from intersection of sphere and plane at a height in between z and z+dz, 
will vary from r+dr to r respectively. Since same size circular profiles will also be obtained in 
between -z and -z+dz, the probability of radius of circular profiles lying between r and r+dr will be 
given by 

PM) = — ‘ 239 > 

R, 

From equations (2.37) and (2.39): 

P(r> = 1W7 am 

The probability density of radii of circular profiles can be written as 

The mean value of radii of circular profiles is given by 

A 1 

= f r. p(r) dr = — 71 R, 

0 4 ' 


r 


(2.42) 



'in .2, Derivation of Sphere size Distribution from a Measured Distribution of 

Profile Sizes 

I 

I 

i 

I 

Suppose that the spheres are distributed inside a cube and the size of the spheres vary from 0 
to Rsmax- Let the probability distribution of radius R have density N v (Rs), so that Nv(R) dR is the 
number of spheres per unit volume whose radii lie in the range (R /R+dR). The probability 
distribution of profile radii have density N A (r) so that Na ( r) dr is the number of circles per unit 
area, whose radii lie in the range (r/ r+dr). The relationship between N A (r) and N\<R) is given by 


N.,(r) dr = Nv( R s ).P Rs P, rt dr dR , 


From equation (2.3 1): 


N, = Nv . 2 .R, 


where Rs is the mean radius of all spheres. 
Hence from equations (2.431 and (2.44): 


(2.43) 


(2.44) 


— T EUkL 
Rs r V^T 7 


(2.45) 


where f(r) = N A (r)/N A is the probability density of profiles and F(R S ) = NvtR)/Nv is the probability 
density of spheres. 


Equation (2.45) is an Abel's Integral equation whose solution is given by 


F(RJ 


d(.j(r)/r) 

4 Rs Rs V dr T _ 

n r, ( r 2 - R; 


(2.46) 


f(r) is the experimentally measurable distribution from the microstnicture obtained on the plane 
of polish. Equation (2.46) can then be used to unfold the distribution of spherical particles dispersed in 
the bulk of a material. However, the above analysis is limited to the spherical shape of particles. For 
non-spherical shapes the solutions become extremely complex and often there is no solution. 


CHAPTER 3 


APPROACH 


Interfaces between particles or grains when cut by the plane of polish yield line segments in 2- 
D microstructure as shown in figure 3. 1 . Transformation from 3-D to 2-D results in the loss of some 
information (typically shape and size distribution) . In some special cases such as spherical particles, a 
complete description of the 3-D structure can be obtained from an analysis of 2-D microstructure, as 
already discussed in section (2.7). However, in most cases this is not possible without making any 
assumptions regarding shape or size distribution of microstructural components. 

In the subsequent chapters, an analysis of the 2-D microstructure obtained from a 3-D structure 
consisting of circular and elliptical interfaces has been done. These interfaces can also be regarded as 
circular or elliptical discs of zero or near zero thickness (thin discs). All such interfaces (or 2-D 
particles) will simply constitute as lines in the 2-D microstructure. The interlaces are assumed to be 
randomly distributed in 3-D space (this assumption is generally true for structures in most real 
materials). The line length distribution is the only experimentally measurable parameter. Using 
principles of stereology. basic stereological parameters describing the 3-D spatial structure have been 
developed in the present work. 

A rigorous analysis of basic stereological parameters has been carried out in order to develop 
the relationships between various 2-D and 3-D parameters,' Stereological analyses of discs 

of circular and elliptical shapes of constant and varying sizes have also been carried out. Practical 
methods have been proposed to extract information about the 3-D structure of these discs from 2-D 
measurements. 

In order to test the validity of various relationships, a distribution of circular as well as elliptical 
discs were generated in 3-D space. This structure is then sectioned to obtain the 2-D microstructures. 
An example of such a structure is shown in figure (3.1). The figure shows that for the case of spatial 






distribution of circular discs of various size, the 2-D microstructure consists of a distribution of line 
lengths. Detailed methodology of such simulations has been discussed in chapter 5. 

From the above mentioned simulated microstructures the following ' ^neasurements\vere made: 

1 . Number of intercepts per unit length , N L 

2. Total length of lines per unit area in the microstructure, L A 

3 . Number of lines per unit area, N A 

4. Line length distribution and computation of statistical parameters, such as, average line length 
and standard deviation of line lengths. 

The above ^measurement^ 1 on the 2-D section are used to obtain some standard 3-D 
stereological parameters such as, surface area per unit volume, Sv and number of discs per unit 
volume, N v . The calculated values can then be compared to the actual values used in the simulation. 
These comparisons were used as a test for the validity' of the simulations (derived in chapter 5 ). 
Correlation between the line length distribution to the disc shape and size distributions have also been 
developed in the subsequent chapters. 



CHAPTER 4 

ANALYSIS OF DISCS 


As stated earlier, when the discs are intersected by the plane of polish, lines of 
different lengths are obtained, as shown in figure 3.1(b). In this chapter an analysis of the 
distribution of the line lengths has been developed to obtain size, size distribution and 
shape of discs in the o-D spatial. The present study is restricted only to circular and 
elliptical discs. For both types of discs several stereological relationships have been 
established. Practical methods have also been developed for the analysis of the 
experimentally measurable parameters of microstructures. 

* 

4.1 Circular Discs 

All the discs are assumed to have a circular shape. The expressions for various 
parameters (tangent diameter, probability of cutting etc.) are obtained from a circular disc. 
The expressions are used to derive the relationship between the distribution of line lengths 
and the distribution of circular discs. First, discs of constant size have been analyzed and 
then the analysis has been extended to discs of varying sizes. 

4.1.1 Discs of Constant Size 

Let the radius of the constant size discs is 'R'. The discs are randomly oriented in 3-D 
space. Figure 4.1 shows a circular disc whose orientation is defined by the angles 9 and o. 

4.1.1.1 Average tangent diameter 


In order to determine the average tangent diameter, the tangent planes are kept 
parallel to the XY plane and the disc is permitted to take all possible orientations, as shown 
in figure 41. From this figure it is dear that there is no effect of angle 6 on the value of 
tangent diameter H(0,<j)) because changing angle fy simply results in the rotation of disc f 
around the Z-axis. H(8,6) can be expressed by the following relation : 


H(Q ,{ JO = 2R\cos(Q -90)\ = 2R\sin(Q )\ (4.1) 


Averaging over all orientations, the average tangent diameter can be written as: (using 
equation (2.23) ) 



(4.2) 


K 


4.1.1.2 Probability of cutting a disc by the plane of polish 


Consider the discs to be contained in a cubic volume, as shown in figure 3.1(a). From 
equation (2.28), probability of cutting a randomly oriented circular disc of radius 'R' by a 
horizontal plane of polish within the cube of dimension 'L' is given by 

P.=y (4.3) 


Substituting H from equation (4.2) in the above expression gives the probability of cutting 
as : 



(4.4) 


4.1.1.3 Average length of lines on the 2-D section 

On the plane of polish, lines of different lengths are obtained as shown in figure 3.1. 
Figure 4.2 shows that the line "mn" is obtained from the intersection of a horizontal plane at 
a height z from the origin and a circular disc of radius R and orientation (9,o). O is the 
origin of XYZ coordinate system and also the centre of the circular disc, 'op' is the 




perpendicular drawn from the origin to the line mn. The perpendicular (op) bisects the line 
run. Length of the perpendicular op is given by. 


Z 

op = 

. 57/2 0 


The length of line mp can be written as 


mp = 



(4.5) 


(4.6) 


Length of line mn is given by j 

K 6.40 = 2 (4.7) | 

V suiv 

i 

] 

It may be noted that length 1(0, (j)) is independent of angle 6. Now the average intersection 
length is given by (using equation (2.4) ).: 

l = -L~j.]l(Q,$).dQ.d<b (4.8) ; 

471 o o I 


kR 

2 


(4.9) 


The average value of reciprocals of line length (harmonic mean) for a circular disc of size R 
has also been given (without derivation) by [39]. 



(4.9 b) 


4.1.1.4 Probability density of line lengths 


Rearranging equation (4.7) as: 



(4.10) 


if e,<j>) = -r 

sinQ 

Differentiating the above equation with respect to z. 

dl_ _ -2z 1 

dz sinQ V/T’sz/fQ - z 2 


(4.11) 


, l.sin' 6 „ 
dz = dl 

4z 


(4-12) 


Negative sign shows that the value of dz decreases with an increases in the value of dl. 
Taking the absolute value of dz 


dz = — dl 


4z 


(4-13) 


Rewriting equation (4.10) as: 


Z 


sinQ 


1 


R- 


4 


(4.14) 


Substituting z from above equation in equation (4.13) . 


dz 


l sinQ 



=dl 

V 

4 


( 4 . 15 ) ; 

! 


Fig (4.3) shows a disc of size ’R' and orientation (0,6). The centre of the disc lies at the 
origin of XYZ coordinate system. A plane perpendicular to Z-axis intersect the disc at a 
height z from the origin. If a plane perpendicular to Z-axis lies within the range of tangent 
diameter (from -Rsin8 to Rsin0), an intersection of the disc with the plane occurs. 




Figure 4.3 : Probability of intersecting of a circular disc. 


Probability of finding a random plane perpendicular to Z-axis, which intersects the disc 
between height z to z+dz from origin . is given by 


P(Z) = 


dz 


Tangent Diameter 


(4.16) 


or 


P(z) 


dz 


2R sin 9 


(4.17) 


The length of lines obtained from the intersection of disc and plane at a height in 
between z and z+dz are 1+dl and 1 respectively. The length of lines are also found between 1 and 
1+dl for a plane intersecting the disc between height -z and -z+dz. So the probability of 
line lengths lying between 1 to 1+dl is given by. 


P(l) = 


dz 

R sinQ 


(4.18) 


Substituting the value of dz from equation ( ), expression for probability' density' of 
line lengths can be written as 

p(l)= =L== (4.19) 

1 - V 

Average length of lines can also be determined by using probability density' of line lengths 



]lp(l)dl 


a 


(4.20) 


From equation (4.19), the equation for average length can be written as: 



(4.21) 



The average length is same as already derived in section (4. 1 . 1 .3). 
Standard deviation (cr) of the line length distribution is given by 


2R 


(1 - l if l 


- J 

1 7 ~ 



= 

r 


dl 


dl - LT 


\R 2 


r 


(4.22) 


4.1.2 Circular Discs of Various Sizes 

The interfaces in any sample would generally have different sizes. These interfaces can 
be approximated by thin circular discs of various sizes. The parameters such as average length 
and standard deviation depend on the size distribution of circular discs. The profile of lines 
obtained from intersection of a random plane with discs is analyzed to obtain the distribution of 
discs. 

4.1.2.1 Average length 

If the size of discs varies from R mn to , the average length of lines can be written 


as: 

R 

'f Probability. density. of .discs * Probability. of .cutting * Average. length. dR 
_ — (4.23) 

R*va. 

J Probability. of .cutting. dR 



Standard deviation of lines can also be given by 


CT = J (l -l )* probability. density. of . lines dl (4.24) 

Rmum 


Example 


For random uniform distribution of discs where the size of the discs varies from 0 to 
Rmax- The value of average length is given by 


R v* 1 4R tcR 
j — — dR 

0 Rmax ff 2 

f — dR 
o nL 


(4.25) 



(4.26) 


4.1.2.2 Derivation of disc size distribution from a measured distribution 

of line lengths 

Size distribution of discs is obtained from the measurements of line lengths. The 
elliptical discs are contained in a cubic volume of unity. Two types of methods are used to 
derive the relationship between n A and n v . 


Method 1. Continuous Method 



Let the probability distribution of the disc sizes have density n%<R), so that n v (R) dR, 
is the number of discs per unit volume whose radii lie in the range (R/ R+dR). Probability 
density of lines lying in the range (1/ 1+dl) after cutting a disc of size R’ is given from equation ( 
4.19). If the probability distribution of the lines have density n A (l), so that n A (l) dl. is the 
number of lines/area, whose lengths lie in the range (1/ 1+dl). The relationship between n A (l) and 
nv(R) can be written by using equation (2.43) 


n A (l)dl 


J m(R) 


4R 

K 


l 

4 r4¥774 


dl dR 


n A (l) 


or, 

-] l n d 3 1 dR 

* Jr 1 -? / 4 


(4.27) 


(4.2S) 


Substituting 1/2 by x in the above equation 


nAD 


2x~. nv (R) 
n i 4r 2 - X i 


(4.29) 


which is an intearal eauation of Abel's type whose solution is given by' i-J 


4 K * 2x r 
nAD = ■' 

K 


R d[nJRUR)dR 

dR 

(r 2 -x : r 


(4.30) 


(4.31) 


riv (l) = 8x J 




(R : -x 2 ) 


.2 J '2 


Integrating the above equation (Integration by parts). 


jiv(l) = 8x{ 


rum T 

m-x 2 r* 

7 

w-x 2 f 


7 ru(R)dR 
Ir(R 2 - x *) ,/2 

dR} 


(4.32) 


or 


?iv(R) = 5x/ 




(D 


(t-x l f 2 * 

fn A (l) 


+ J 




;/2 


7 n.A 

» /r / 2 - r r 

dlj 


(4.33) 


For solution of practical problems replace °° by W x /2, where is the maximum line length 
obtained in the micro structure and also replace x by 1/2. If the length distribution iu(1)i is 
known, the above equation can be used to unfold the distribution of circular discs ny(r) 


Examples 

1 . The following equations can be derived for random distribution of discs, where size of 
discs varies between 0 to Rm* 



Say, n v (r) = C = Constant 


Ua (1) = 



71 1/2 


nv(R)dR 



1C 

n.\(l) - — [ log( seed + tanQ ) ] 

K X =!/2 


1C 2R \4R - *=*- 

nM = — [ ] 

71 l V / *« i/a 


m 


■(R) = 


n n A (l) 


l[log(2 


Rn 


4R z n 


-D] 


2. If the distribution of line lengths is given by equation: n A (I) = 1 e' such that : n A (0) 
and n A (=<>) = 0 . Substituting the value of n A (I) in the equation ( 4.33) 


riv(R) = 8x[ \ 


e'dl 7 le‘dl 7 e dl 


r * t: I* . e -t 

J / 7 2 2 J / ;2 2 \I /2 / - \ i/2 

*( l ~ X ) *(l -X ) M / ' X ) 


Thus .the size distribution of the discs is given by 




In general, equation (4.33) will need to be solved numerically. Often the approach of 
assuming continuous distributions can be tedious and therefore, it may be easier to discretise 
the distribution, as discussed below. 

Method 2. Discrete Method 

The actual distribution of discs is approximated as discrete distribution where the disc 
sizes are divided in se\ eral classes. All the discs in a particular class can then be assumed to be 
the same size. Let the size of the discs varies from 0 to Rmax- The discs are divided in 'M' 
classes and the width of each class is 

AR = R™x/M (4.34) 

So the i lh class will contain all the discs whose radii lie in between (i-l)AR and iAR and the 
value of i varies from 1 to M. The maximum length of line which can be obtained from the 
distribution of the discs will be the diameter of the largest size disc. So the length of lines will 
vary between 0 and ZR^. 

The length of lines are also divided in M classes and the width of each class is 

2AR = 2 Rna/M 


So the j* class will contain all the lines whose length lie in between (j-l)2AP and j.2AR.and the 
value of j also varies from 1 to M. 



Assumption 

The distribution of disc sizes is assumed to be discrete i.e. a particular class has a 
particular size only. The size of the discs in a particular class has been chosen as the largest 
size in that class. However, the intersection lengths are considered to assume a continuous 
distribution which are also then broken up into classes, as given in table 4.1: 


Table (4.1) Distribution of discs and lines in various classes 


Class No. 

Size 

Numerical Density 

CIRCULAR DISCS 

1 

AR 

Nvi 

2 

2AR 

n v2 

3 

3AR 

n v3 

i 

iAR 

N^ 

m 

mAR 

Nym 

Profile (Lines) 

1 

Oto 2AR 

Nai 

2 

2ARto 4AR 

Na2 

3 

3 AR to 6AR 

Nas 

j 

2(j-l)ARto 2jAR 

NAi 

m 

2(m-l)AR to 2mAR 

Nam 


Since the size of the discs in the i* class is iAR,the probability of cutting a disc in the i class is 
given by (from equation (2.28) ): 




Pi = 


4iAR 

tiL 


( 4 .;- 6 ) 


ForL= 1, 


Pi = 


4iAR 


iz 


(4.37) 


Probability of lines lying in the j* class after cutting the discs of ith class (i > j) is given by 
(from equation (4. 19) ) 


d . _ Idl 

4R4WT74 


(4.38) 


IN the above equation the value of 1 is taken as the mid point of a particular class and therefore, 
Pi J becomes : 


pi = 


(j-l)2 AR + J 2AR 

_ 2 

I ~ (2j - 7 Jar ~ 

4iAR , (iAR r-[~d — ]- 

V 4 


(4.39) 


or. 


Pi 




(2j-l ) 


i-\4C-[2j-l ■ f 


(4.40) 


The number per unit area of lines in j* class contributed from the discs of ith class is given by 


P, * Nv (i) * Pi 


(4.41) 


Therefore, the total number of lines per unit area in j -2 class is given by 



3 


n Aj = P m * Nv „ * PL + P m _y * N v „., * PL; 

. (4.42) 

+ + Pj * N Vj * />', 

or, 

v m 

Ha =1P ; N v P/ (4.43) 

i=j 1 

Above equation can be expanded to ‘m’ simultaneous equations with ‘nf unknowns, i.e, n V i. 

Hv 2 5 - * - - . Uvm '• 

Forj = m 


N a . = 


4M 


7C 


M * i\r v 


2/71 - 1 


m 


44rn-l 


or 


7V„ 


7c44m-l 


A R_Nr w 


(4.44) 


(4.45) 


Forj = m-1 





4 m 


71 


AjR * N„ 


* , 2m ~l + — V AR * .V 

v3-v4/h-3 7t 


+ (2m - 3j 

(m- 1 )*j4m — 5 


N a . 


4(2m - 3) 


-J3ii-J4m - 3 


AR N v 


4(2m -3) 


K -44m - 5 


■AR Nv... (4-47) 


The number of lines in other classes can also be written in similar fashion. The set of m 
equations given by equation ( 4.43) can be easily solved to determine m unknowns. Thus 
knowing the values of N A1 - N^, one can calculate the values of Nvi ...... N Vm . 


4.2 Elliptical Discs 


All the discs are assumed to have elliptical shape. The stereological relationships for 
various parameters are obtained in the same fashion as were obtained for circular discs. First 
elliptical discs of constant size are considered. Subsequently discs of varying sizes are analyzed. 
However, the shape of elliptical discs are assumed to be a constant, i.e., the ratio of major axis 
to minor axis is kept constant. 

4.2.1 Discs of Constant Size 

The distribution of discs is assumed to have all the discs of single size. A disc from the 
distribution is chosen randomly. Let he semi-major and semi-minor axes of the elliptical disc are 
'a' and 'b' respectively. The orientation of the disc is given by 9,q and c, as shown in figure 4.4. 

It is necessary to define an extra parameter (q) for the orientation of elliptical disc because the 
rotation of disc about its normal changes the position of the disc. While the rotation of a 
circular disc about its normal has no effect on the position of the disc. The following 
parameters are obtained for the disc mentioned above. 

4.2.1. 1 Average tangent diameter 

Tangent Diameter of the elliptical disc is taken as the distance between tangent planes 
parallel to XY plane, as shown in Figure (4.4 ). There is no effect of angle o on the tangent 
diameter H(0.b.q) because change in <p merely results in the rotation of discs about the Z-axis. 
The tangent diameter of an elliptical disc having orientation 0 , 0 and c can be written as 

H(Q, <j>, £) = H(Q = j, CJ sin 0 (4-55) 

where. H(0 = m/2, q) is the maximum value of tangent diameter at 9 = rc/2 and at a constant value of 

5 




In order to determine the value of H(6=n/2 , q), the ellipse is fixed in XYZ- 
coordinates in such a way that its major and minor axes lie along X and Y axes respectively. 
The equation of ellipse in XYZ-coordinates is then given by 



and 
Z~ 0 


( 4 . 56 ) 


Instead of rotating ellipse by (0=7t/2, !;) in anti-clockwise direction, the axes are rotated 
clockwise by the same amount. After the rotation, XYZ coordinate system is changed to 
X'Y'Z'. The parallel planes which are tangent to the disc and perpendicular to the new Z axis 
represent lines AB and AB' in XY-plane as shown in figure ( 4.5). 

Equation of line AB in XY coordinates is given by 


y ~ tan (90 + t,)x + c 

or 

y = -cot t,.x + c 


(4.57) 


where, c is the intercept of the line AB on the Y-axis. 

Substituting equation (4.57) in equation (4.56) and rearranging : 

x 2 (b : + a 2 cot 2 t) + x (-2 a 2 c.cot q) + a 2 (c -b' ) = 0 (4.58) 


Since the line AB just touches the ellipse, the solution of equation (4.58) should give a 
single value of x and y. The two roots of the equation (4.58) can be made equal by using the 
following property 

4aV cot 2 £, = 4(b 2 + a 2 cot 2 !;). a 2 (c 2 -b 2 ) 


or 


c = (a 2 cos^s + b 2 sin'Q f~ / sine 


(4.59) 



Substituting c from above equation in equation (4.57), the equation of line AB becomes : 


y - -cot^.x + (a cos\ + b 2 sin %)' u / siivt (4.60) 


Figure ( 4.5) shows that the line CD, which is perpendicular to the pair of lines AB and 
AB' and passing through origin, is along the new Z-axis. This line intersects the lines AB and 
AB' at C and D respectively. The tangent diameter of the ellipse having orientation (8 = tc/ 2, 
q) is nothing but the length of line CD. The equation of the line CD is obviously given by 


y = tan 


(4.61) 


From equations (4.57) and (4.59), the coordinates (x,y) of point C is given by 

.x = - cos ^ ( a cos £ + b 2 sin C ) U ~ (4.62) 

and 

y = sin £,( a 2 cos + b~ sin C f‘ (4.63) 

Length of line CD. i.e.. H(0=7t/2 , Q, is simply given by 

H(Q = 7C / 2,\) =2 ( x 2 + v' ) W1 —( a~ cos' +-b~ sin £, ) (4.64) 

Thus, from equation (4.55), tangent Diameter H(9,cj),q) is given by 

H(Q,fo&) = 2(acos 2 £, + b' sin C ) .sinQ (4.65) 

The average tangent diameter ( H ) can be written as: 


CENTRAL UBRARf 



H = 


[ J ]H(Q^,^)dQd^ 


J J J dQdtyd^ 


0 0 0 


* '7C 

J 1 ( a cos' ^ + b 2 sin ' dq sinQ dQ 


(4.66) 


Suppose 


2 - j (a co s' 5 + b 3 sin 2 5 f* 

0 

= (~ L r1{l + ~ 7 cos2^/ ,1 d^ 

2 o a+b 2 


(4.67) 


The solution of the integral® can only be written in the series form[ 20 ]. 


/ = 2k(SL±£ /'7i..iL..i£ 

2 16 7“ 


t 4 -■ 


where t = (a 2 -b : )/(a 2 +b 2 ) 

Therefore, the average tangent diameter is given by. the following relationship : 



K 



\ l /2 / 7 ^ "* ^2 

) t~-—r 

? 9* J 




(4.68) 


(4.69) 


For the case when a = b = R (i.e., circular disc), the equation of average tangent diameter for 
elliptical disc reduces to equation (4.2) which is identical to the equation of average tangent 
diameter derived for circular disc of radius R. 


4.2.1.2 Probability of cutting a disc by the plane of polish 



From equation (2.28), probability of cutting a randomly oriented elliptical disc by a horizontal 
plane of polish within a cube of dimension L is given by 


*=± ( l±]L^n-L. 15, 


L Til/ 2 


ra-~r- 

2 2 




(4.70) 


where t = (a 2 -b 2 )/(a 2 +b 2 ) 


4.2. 1.3 Average length of lines 


The elliptical disc whose semi-major and semi-minor axes are 'a’ and 'b' and has the 
orientation (0, <(>,£) with respect to X-Y-Z frame, is cut by a plane perpendicular to Z-axis at a 
height z from the origin as shown in figure (4.6). 

Consider the frame X'Y'Z', which is obtained from XYZ frame by following operations 
in sequence. 

1 . Rotation of Z-axis by an angle 0 in anticlockwise direction 

2. Rotation of new Y-axis by an angle <j) in anticlockwise direction 

3. Rotation of new Z-axis by an angle q in anticlockwise direction. 

Since the orientation of the normal of the disc is defined by (0,o,q) in XYZ frame and the 
orientation of X’Y’Z’ can also be given by (0,<j>,q) with respect to XYZ frame. So, the normal of 
the disc v/ill lie along Z’-axis and the disc will lie in X’Y’ piane. as shown in figure 4.7. The 
equation of ellipse in X'Y'Z’ frame is given by 



a b 


In the figure (4.7 ), axis Z' is coming out of the plane of paper . Line AB is the intersecting 
line between elliptical disc and sectioning plane and line OD is perpendicular to line AB. The 
length of line OD is given by 



z 


(4.72) 


OD = - 

sin 9 

Equation of line AB in X' Y' frame is given by 

y' = tan(90 + ^)x' + — - (4.73) 

sinQ sinq 

Substituting equation (4.73) in equation (4.71) : 

( m a + b 2 )x f 2 + ( 2m a 2 c)x ' + a(c-b 2 ) = 0 (4.74) 

where m = tan (90+q) and c = z/sin0 sine. 

Let xi and X 2 be the two roots of the quadratic equation (4.74), one can write the following 
identity : 


( Xi - X 2 )‘ - ( Xi + X 2 )‘ - 4 xi Xi (4.75) 

Using the properties of summation and multiplication of roots® .one can write the following 
equation by using equation (4.75) : 


(x ,-x 2 f = 4 a 2 b 2 (a 2 m 2 + b 2 -c 2 )/(m : a z + b 2 ) 2 
From equation (4.73), one can write the following equation 



m m 

where m = tan (90+q) and c = z/sin0 sinq. 


(4.76) 


(4.77) 


A quadratic equation in terms of y ! can be obtained from equations (4.71) and (4.77) : 


If xi and x; are die two roots' of equation air + bx + c = 0, the following identities hold good 
xi -r xz = -b/a; 
and 

X1.X2 = c/a . 


( m a +b)y' 2 - (2 b 2 c)y' + b 2 (c~a) = 0 (4.78) 

If yl and y2 are the roots of the above equation .then the square of the difference of the roots is 


given by 


( y i - y 2 f 


( m a + b 2 f 


(4.79) 


Therefore, the intersection length between the ellipse and the plane of polish , i.e., the length AB in 


figure 4.7 is given by 


WA&) = {(x,-xj + (y r yj } 1 ' 3 


(4.80) 


>§>£>) - : , 7Z ; jr--{a*sin 2 Qcos 2 £ ) + b~smQsm£ ) -z'} l/ ‘(,<>]\ 

sinQfacosq + b'sinq) (4541) 

By putting a = b = R (i.e., case of circular disc ) equation (4.81) reduces to the equation (4. 10) 
which was derived for circular disc. 

Average length of lines in the microstructure is obtained by averaging over 9.o and q as follows 


nr 2% HfQ.o.c, ) 

j j j IdzcTqdQ 

0 0 -H(Q. o£) 

7t 2zH(B.O,c > ) 

j J J dzdidB 

0 0 -HtQ.o.c,} 


(4.82) 


7 is me, 04 ) 

1= ff j Idzd^dQ 

4n : H(Q,§£)a 0 -mo.oAi 


(4.83) 


4.2. 1.4 Probability Density of Line Lengths 



Rearranging equation (4.81) as : 


z = 


sin 9 
2ab 


{ a cos 2 £> + b 2 sin £ f'~ { 4 a b 2 -f( a 


2“ , ,2 . 2p jlf2 

cos 's + b sm s / 


(4.84) 


Differentiating the above equation with respect to 1 


d - + sinQdl 

Z ~ 2ab{4 a 2 b 2 -l 2 (a 2 cos% + b 2 sin%)} ,n 

The following expression of probability of lines lying between 1 to 1+dl can be obtained in an 
analogous manner as in the case of circular disc (See section 4. 1.1.4). 

P(l) = 2dz/H( 9,<j>^) (4.86) 


or 


P(WV = 


(_ a_ cop sin l )ldl 

2ab{4 a b 2 - 1 2 ( a cos 2 £> + b 2 sin sin 


(4.87) 


where P(l){q} is the probability density function of line lengths at a particular value of q. The 
probability density of line lengths can be obtained by integrating over the range of % is given by 

P(l) = -Lj P(l)f^) ^ (4.88) 

271 o 

4.2.2 Elliptical Discs of Various Sizes 

The interfaces in any sample generally have different sizes. In some cases these 
microstructures are found to have a shape close to the elliptical shape. In order to analyze the 
distribution of such microstructural elements, the microstructure can be approximated by thin .j 
elliptical discs of various sizes. The ratio of semi-minor and semi-major axes (b/a=k) is assumed to I 
be constant (i.e, the shape of elliptical discs is kept constant). The discs are analyzed using a similar 
approach as used for circular discs. 



4.2.2. 1 Deri\ ation of disc size distribution from a measured dis tribution of 

profile sizes 


When a cubic volume filled with elliptical discs of various sizes is intersected by a random 
plane, the lines of different lengths are obtained on the plane of polish. Distribution of the discs inside 
the cubic volume is obtained from the measurements of line lengths distribution. The elliptical discs are 
contained in a cubic volume of unity. Two approaches are used as discussed below' : 


Method 1. Continuous Method 


Let the probability distribution of the semi-major axis have density n v (a), so that ny(a) da, is 
the number of elliptical discs per unit volume whose semi major axes lie in the range (a/a+da). From 
equation (4.69) and noting that cube dimension is unit ( L = 1), the probability of cutting the elliptical 
disc of semi-major and semi-minor axes a and b respectively, is given by 


P = 


71 


2 J. 

■ra- 


4a ,1+ k‘ an r , t 


15 


16 2 


10 




(4.89) 


where k = b/a. 

Probability density of lines King in the range (1/ 1+dl) after cutting a disc of size a is given by equation 
(4.88). 

p(d = 2-7 ( 4 . 90 ) 

2n o 

Let the probability distribution of lines have density n A (l) per unit area, so that n A (l) dl is the number of 
lines/area, whose lengths lie in the range (!/ 1+dl). The relationship between n A (l) and n v (a) can be 
written by using equation ( 2.43) 

n A (l)dl = ] n ,(a).H.P(l)da (4.91) 


or 



} 


tia(I) — —~jn v (ci) 


2tz 1/2 


*?- ( 1+ k f2 n r lo , 

_ ) -M- — — -t - 

71 2 16 ? J0 


* j' [ a_ co£_ £, + b 2 sin_ £, J/ _ 

o 4 a 2 b 2 -l 2 ( a cos 2 £ + sin 2 £ j / /J 


(4.92) 


The above equation can only be solved numerically by assuming some kind of a distribution for N v (a). 
As will be shown below, by assuming discrete distributions instead of the continuous distributions, one 
can obtain a much simpler solution for practical problems. 

Method 2. Discrete Method 

The actual distribution of discs is approximated as discrete distribution where the discs are 
divided in several classes. All the discs in a particular class are assumed to have the same size. Let the 
size of the disc varies from 0 to a^x. And the discs are divided in M classes such that the width of each 
class is 


la = — (4.93) 

M 

Therefore, the i lh class will contain all the discs whose size lie in between ii-l).Aa and i.Aa and the value 
of i varies from 1 to M. 

The length distribution obtained from these discs will have the length of lines varying between 
0 and 2Aa max . The maximum length of line being the length of major axis of the largest size disc. The 
length distribution of lines are also divided in M classes such that the width of each class is 

2ia = (4-941 

M 

So the j* class will contain all the lines whose length lie in between 2(j-l )-Aa and 2j.Aa_ and the value of 
j will vary from 1 to M. 



Assumption: 


The distribution of the disc sizes is assumed to be discrete i.e. a partic ular class has a particular 
size only. The size of the discs in a particular class has been chosen as the largest size in that class. 
However, the intersection lengths are considered to assume a continuous distribution which are also 
then broken up into classes, as given in table 4.2. 


Table 4.2 . Distribution of elliptical discs and lines in various classes 


Class No. 

Size 

Numerical Density 

ELLIPTICAL DISCS 

1 

Aa 

Nvi 

2 

2Aa 

n v2 

3 

3Aa 

Nv3 

i 

iAa 


m 

mAa 

Nvm 

Profile (Lines) 

1 

0to2Aa 

Nai 

2 

2Aa to 4Aa 

Naz 

3 

3 Aa to 6Aa 

Nab 

j 

2(j-l)Aa to 2jAa 

Nai 

m 

2(m-l)Aato2mAa 

Nahs 


Since the size of the discs in the i 8 * class is iAa, the probability of cutting a disc in the i class is given by 
(from equation (4.70) ) : 


4.i,Aa ,1 + K 2 \m n l_ a > 

jt ' 2 J i 'i6~r 


(4.95) 


Probability of lines lying in the jth class after cutting the discs of ith class (i > j) is given from 
equations (4.87) and (4.88) 




pi - ~ ll 2 J f l'05'C ^ K 2 sin ^ 

2k.LK o {4i 2 K 2 - (2j -lf( cos\ + K 2 sm 2 l)}' 

The number of lines in j* class contributed from the discs of ith class is given by 

Pi * NM) * Pi 

Total number of lines per uniy area in j* class is: 


(4.96) 


(4.97) 


Naj = Pm * N Vm * Pm + P m -1 * iVvw * Pil + + Pj * Nr, * P j (4.98) 


or, 

= X Pi * Nr, * ?/ (4-99) 

<=/ 

For different values of j (from 1 to M) equation (4.99) will give m simultaneous equations, 
which can be solved for M unknowns ( N VI , N Y2 , N Vm )• 




CHAPTER 5 


MICROSTRUCTURAL SIMULATION 

Simulation was used for evolving a model which represented the real sample. Since 
it is difficult to perform the experiments on real samples, microstructures containing 
circular and elliptical discs were simulated. The simulated microstructure were used for 
testing the equations developed in the previous chapter (see appendices C-l to C-4). 

Consider a cube of dimension 1, containing circular or elliptical discs. If the cube has 
hlv number of discs per unit volume, then Wv is the number of discs inside the cube of 
dimension 1. In order to simulate a 3-D structure nv number of discs are then randomly 
distributed (which also have random orientations ) inside the unit cubic volume. Finally, 
random sections of the cube are taken to generate microstructures on the plane of polish. 

A disc of size ai (semi-major axis) with orientation (9i, bi, gi) is generated at any 
random position inside the cube. Another discs of size a 2 (semi-major axis) with orientation 
(02, § 2 , C 2 ) is also generated at some other location inside the cube. Suppose the centres 
of the first and second discs are at (xi, yi, zi) and (x;, y 2 , zi) respectively. 

5.1 Generation of Distribution of Discs in 3-D Space 

j 

■ j 

In order to generate a random position and random orientation of a disc the 
following procedure was adopted : J 

• A random position of the disc was obtained by generating o uniform random numbers 
corresponding to the x,y,z coordinates which denoted the centre of the disc. 



. For random orientation 0 and $ were generated for circular discs and 0,<j) and £ were 
generated for elliptical discs. It should be noted that $ and $ are uniformly distributed 
from 0 to 2k and in the case of 9 , cosQ is uniformly distributed over -1 to +1 (i.e. ; 9 going 
from 0 to 7t ). 

During the generation process described above it was ensured that no overlap of discs 
occured. For this the conditions of overlap between two discs have been discussed in the 
next section If an overlap occurs, the position and orientation of the last disc was generated 
again. 


5.1.1 Condition for Overlap of Two Discs 

Let the positions of the centres of two discs be (xi,yi,zi) and (x^y 2/ z 2 ) respectively. 
The orientation of the two discs are given by (9i ,<jh) and (62 ,02) for circular discs and 
(0i,<t>i, £,i) and ( 0 2 ,62,02) for elliptical discs. For the case of circular discs, the radii of the two 
discs are denoted by Ri and R2 respectively. For the case of elliptical discs the semi-major 
and semi-minor axes of the discs are denoted by (ai,bi) and (a2,b 2 ). 

Suppose 

Ax = x 2 - Xi 
Ay = y 2 - yi 
Az = z 2 - Zi 

Condition 1. 

If | Ax | or | Ay | or | Az I > ai + a 2 , there is no overlap between the discs. 

Condition 2. 

If di <= bi and d 2 <= b 2 the two discs overlap each other. Otherwise , there is no 
overlap between the discs. It may be noted that for the circular discs, ar = bi = Ri and a 2 = b 2 
= R 2 . The expressions for dr and d 2 axe complex and therefore they are given in appendix-B. 

Here di and d : are the distances from the centers of first and second disc to the 
intersecting line between the disc planes respectively. 



5.2 Generation of Microstructure on the Plane of Polish 


In order to generate 2-D microstructures exhibiting lines of varying lengths, the cube 
of unit dimension containing randomly distributed discs was sectioned randomly by 
planes parallel to the XY-plane. The position of the plane of polish along the Z-axis is 
chosen randomly. 

The length of line intersected between a disc and the random plane is given by 

i = L(x x - ny ) 2 + ( y\ - y 2 ) 2 ] 1/2 

where (xi,yi) and (x 2 ,y 2 ) are the end points of the line. The expressions for xi , yi , x 2 and y 2 
are derived in appendices A-land A-2. 


For each simulated distribution of discs, the cube was sectioned randomly 10 to 20 
times to obtain 10 ,to 20 microstructures of lines of the type shown in figure 3.1(b). 


CHAPTER 6 
DISCUSSION 


Using the methodology discussed in the previous chapter, discs of various shapes (i.e., different 
ratios of the minor axis (b) to the major axis (a) of ellipse, k=b/a) and sizes were generated in a cubic 
volume. The cube was sectioned by random planes (or the so called plane of polish) and on each 
section a microstructure consisting of straight line segments of varying lengths were obtained. From the 
microstructure consisting of straight lines, the values of 2-D stereological parameters : number of 
intercepts per unit length(N L ), length of lines per unit area (L A ), number of lines per unit area (N A ), 

average line lengths( t ) and length distribution of lines were obtained. Using the values of these 
parameters for the simulated microstructures, 3-D stereological parameters were calculated, namely, 
surface area of discs per unit volume (S v ) and number of discs per unit volume (N v ). Tables (6.1) and 
(6.2) compares the calculated values of Sv and N v with their actual values for structures consisting of 
discs of different size distribution and shapes (circular shape with k=l to elongated ellipse shape with 
k=0.1). 


Table (6.1) : Data for surface area per unit volume 


Disc Shape 

k 

Actual surface area per unit volume 

Sy 

(mnr/mnT) 

Calculated Sv 

Sv = 2*N : 

(mnr/mmf) 

Calculated Sv 

Sv = 4*L a /z 

Discs ot constant size (semi-major axis = 0. Imrr 

i; N v = 400) 


1.0 

1.57 

1.51 

1.49 

f; 0.9 

1.41 

1.38 

1.40 

0.8 

1.26 

1.18 

1.10 

0.7 

1.10 

1.11 

1.00 

0.6 

0.94 

0.93 

0.88 | 

* 

0.5 

0.78 

0.70 

0.72 j 





Table 6.i continued : 


Disc Shape 

Actual surface area per unit volume 

Calculated S\- 

Calculated Sv 

k. 

Sy 

Sv = 2*N l 

S v = 4*L.v'ri 


(mm 2 / mm'") 

(mnr/mnf) 

(ramW) 

0.4 

0.63 

0.60 

0.58 

0.3 

0.47 

0.47 

0.42 

0.2 

0.31 

0.34 

0.32 

0.1 

0.16 

0.18 

0.15 

Uniform random distribution ot disc sizes (semi-major axis between 0 and 0. 1; \~ v = 400 ) 


0.50 

0.49 

0.44 

0.9 

0.49 

0.44 

0.45 

0.8 

0.42 

0.40 

0.38 

0.7 

0.33 

0.31 

0.36 

0.6 

0.31 

0.28 

0.30 

0.5 

0.24 

0.22 

0.22 

0.4 

0.20 

* 0.18 

0.18 

0.3 

0.14 

0.13 

0.16 

0.2 

0.11 

0.13 

0.11 

0.1 

0.06 

0.06 

0.05 j 

Normal distribution of disc sizes (p = 0.2 ; <j = 0.06; N v = 100) 

1.0 

14.26 

13.20 

12.94 

0.9 

12.19 

12.08 

1 1,67 

0.8 . 

11.42 

10.S6 

10.89 

0.7 

9.72 

9.31 

9.33 

0.6 

7.55 

7.54 

7.00 

0.5 

6.91 

7.02 

6.38 

0.4 

6.11 

5.99 

5.88 


4.23 

4.08 

3.92 j 

Muffin 

2.81 

2.77 

j 2.67 

0.1 

1.29 

1.31 

! 1.22 t 

; | 
















Table 6.2 : Data for number of lines 


per unit area 


Disc shape 

k = b/a 

Actual number of discs per unit volume 

N v 

Calculated N v 

n v =nv7/ 

Discs of constant si 

ze (semi-major axis = 0.1mm) 

1.0 

400 

425 

0.9 

400 

412 

0.8 

400 

392 

0.7 

400 

409 

0.6 

400 

410 

0.5 

400 

405 

0.4 

400 

391 

0.3 

400 

355 

0.2 

400 

432 

0.1 

400 

395 j 

* 

Uniform random dis 

:ribution of disc sizes (semi-major axis between 0 and 0.1) 

1.0 

400 

j 391 

0.9 

400 

415 

0.8 

400 

420 

1 

0.7 

400 

400 

0.6 

400 

431 1 

i ! 

0.5 

400 

| 374 

♦ 0.4 

400 

j 383 | 

0.3 

400 

j 386 j 

0.2 

400 

405 i 

0.1 

400 

j 410 j 

Normal distribution of disc sizes (u - 0.2 ; a - 0.06) j 

1.0 

100 

i 105 j 

0.9 

100 

! 101 



















Table 6.2 continued : 


Disc shape 

Actual number of discs per unit volume 

Calculated Nv 

-5 

11 

N v 

N V = N jjff 

0.8 

100 

104 

0.7 

100 ' 

99 

0.6 

100 

95 

0.5 

100 

101 

0.4 

100 

106 

0.3 

100 

103 

0.2 

100 

96 

0.1 

100 

95 


It is clear from tables 6.1 and 6.2 that there is good agreement between the calculated and 
actual values of Sv and N v . Since, equations (2.13) and (2.31) (used for calculating S v and NV) are 
fundamental relationships in stereology, the methodology used for generating simulated 

microstructures is valid. The value of mean tangent diameter. H (used in calculating NV) was 
estimated from equations (4.2) for circular discs (k=l) and equation (4.69) for elliptical discs (k<l). 

From the self-consistent results, it may be concluded that the above equations for H are valld - 

The mean ( / ) and harmonic mean (1) of line lengths can also be obtained (-measured'’) from 

the line segments on 2-D micro structures. Equations (4.9) and (4.9b) were used to calculate the mean 
and the harmonic mean ofline lengths respectively for circular discs. A comparison of the measured 
and calculated values are shown in table (6.3 ) 











Table 6.3 : Data for mean and harmonic mean of line lengths 


Constant size (radius) of 

circular discs 

(mm) 

“Measured” 

from 

microstructure 

1 

(mm) 

Calculated from 

equation (4.9) 

7 = ;uR/2 

(mm) 

Measured” 

from 

microstructure 

(1/mm) 

Calculated 

from 

equation 

(4.9b) 

fl) =^'4R 

(1/mm) 

0.1 

0.1599 



0.1571 

8.1117 

7.8540 

0.2 

0.2918 

0.3142 

4.5637 

3.9270 

0.3 

0.4794 

0.4712 

2.4137 

2.6181 

0.4 

0.6423 

0.6283 

1.6984 

1.9635 

0.5 

0.7720 

0.7854 

1.6370 

1.5708 

0.6 

0.9436 

0.9425 

1.2509 

1.3090 

0.7 

1.1311 

1.0996 

1.1040 

1.1220 

0.8 

1.2383 

1.2567 

0.9552 

0.9817 

0.9 

1.3713 

1.4137 

0.9115 

0.8727 ~| 

5 

1.0 

1.5312 

1.5708 

0.7321 

0.7854 


From the above table, it is clear that the equations (4.9) and (4.9 b) hold good for circular discs 
simulated in 3-D structure within reasonable limits. However, the average value obtained for 
reciprocal lengths (harmonic mean) is more prone to measurement error because small errors in the 
smaller line segments can be magnified many-fold. While the same error will not be magnified in the 
case of mean lengths. If the chances of error in the measurement are higher or the profile has man\ 
small line segments, it is recommended that for the calculation of the average disc size, equation (4.9) 
be used instead of the equation (4.9b). 

The 3-D structure simulated for certain distribution of interfaces (discs) can be unfolded using 
the length distribution of lines obtained from 2-D microstructure by using the methodologies described 













The 3-D structure simulated for certain distribution of interfaces (discs) can be unfolded using 
the length distribution of lines obtained from 2-D microstructure by using the methodologies described 
in chapter 4. Figures (6.1), (6.2) and (6.j) compare the actual distribution of discs as simulated in 3-D 
space with the unfolded distributions of discs. The nature of the unfolded distribution remains more or 
less same as compared to the actual distribution shown in these figures. The values of the calculated u 
and a are found to be very close to the actual values as shown in table 6.4. 

Table 6.4 : Actual and calculated data for u and o 


Disc shape 
k=b/a 

Actual distribution parameters 

from 3-D structure 

Calculated distribution parameters 
after unfolding the micro structure 


u 

a 

u 

a 

Discs of co 

nstant size | 

1.0 

0.1 

0 

0.0.995 j 0.0005 

0.9 

0.1 

0 

0.0912 

0.0103 

0.8 

0.1 

0 

0.0904 | 0.0293 

0.7 

0.1 

0 

0.0913 

0.0408 

0.6 

0.1 

0 

0.0892 

0.0612 

Uniform rai 

[ — — 

idom distribution of disc sizes 

1.0 

1.0 

0.3 

1.0313 | 0.2620 

1 ■ 

0.9 

1.0 

0.3 

0.9316 1 0.3351 

1 

0.8 

1.0 

0.3 

0.9447 | 0.2665 

0.7 

1.0 

0.3 

0.9532 j 0.3116 

Normal di: 

. l 

5 tribution of disc sizes j 

1.0 

1.45 

1.5 

1.4978 , -.j 1.5186 

— : — i 

0.9 

1.45 

1.5 

1.3556 

j 1.3382 

i — 1 

0.8 

1.45 

1.5 

1.4616 j 1.5750 

0.7 

1.45 

1.5 

1.5612 1.6823 



o.s 


0.6 


0.4 






-U 


0.02 0 04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
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Figure 6.1 : Distribution of discs of constant size 

(a) i-Actual distribution for both k=l and k=0. 6 

(b) unfolded distribution for k=l 

(c) unfolded distribution for k=0. 6 
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Figure 6.3 : Normal distribution of disc size 

(a) Actual distribution for both k=l and k=0.6 

(b) unfolded distribution for k=l 

(c) unfolded distribution for k=0.6 




It is not always possible to unfold the distribution as one may not have the complete shape 
information (which is usually lost) from the 2-D microstructure. In order to describe the distribution 
of elliptical discs completely, it is necessary to know the distribution of both semi-major and semi- 
minor axes or the distribution of one of the axis and their ratio. This is a bi-variate problem and in 
order to solve it, bi-variate information is needed, but the section profile provides only uni-variate 
information ( lines are uni-variate ). So, from the profile of lines, it would generally be impossible to 
determine the distribution of both the axes. In order to determine the distribution of both the axes of 
the elliptical discs, some extra information like the nature of the distribution of discs , the disc shape (k) 
etc., is required. Equation (4.99) consists of m equations in terms of m-1 unknowns 

(Nvi,Nv 2 , N Vm and k) which can only be solved for m unknowns. In order to reduce the number 

of unknowns to m, some techniques are proposed to determine the shape factor, k. The methods 
assume a constant (as a first approximation an average value) shape factor, k. as inclusion of a 
distribution of shapes would introduce far too many unknowns and the problem would become 
intractable. 

Equation (4.99) can be solved for different values of the shape fact or, k to unfold the various 
possible size distributions of discs. The different size distributions exhibit different characteristics for 
the appropriate value of k as give in table 6.5. 

Table 6.5 : Characteristics of various distributions 


Distribution of discs 

Characteristics 

1 . Constant size 

.411 the discs lie in the last class 

2. Uniform random 

The discs are evenly distributed among the classes 

3. Normal 

The discs are distributed normally 

4. Lognormal 

Log normal distribution for the discs is obtained | 

5. Others 

Distribution similar to the initial distribution appears for the appropriate j 


value of k j 







Figures 6.4 to 6.7 show various plots for probability density' of lines versus length of lines. A 
similar plot can also be obtained from the 2-D profile of lines. The plot thus obtained can then be 
superimposed on these figures. The value of k for which a plot in the respective figure is found close 
to the plot obtained from the 2-D profile, is the appropriate value of k. However, such a 
superimposition may result in multiple possibilities. 

The coefficient of v anation defined as the ratio, 0 / u. was calculated for the unfolded size 
distributions of discs over a range of k varying from 1.0 (for circular case) down to small values. The 
variation of the calculated ov'u. was plotted as a function of k. as shown in figures (6.8), (6.9) and 
(6.10). Figures (6.8), (6.9) and (6.10) correspond to actual distributions: normal log-normal and 
gamma distributions respectively. In each figure there are several plots corresponding to actual shape 
of discs given by k=1.0, 0.9. 0.8 and 0.7. It may be noted that the curves in these figures are plotted 
only up to a certain minimum value of k. The calculations showed that beyond a certain minimum k , 
the Ny of discs in certain size classes become negative. Since negative values of Ny has no physical 
significance, the cut off value of k can be determined. Thus, it may be concluded thai it is possible to 
determine a range of values of the shape factor, k from the distribution of line lengths. 

The values of Ny in some classes start appearing negative below a certain value of k as 
discussed above. In the figures 6.8 to 6.10. the only values of k for which all the classes have some 
positive value of Ny are shown in appendix D-l. For a value of k larger than the actual \ alue. one 
needs less number of discs in higher classes to get the required number of lines and their contribution to 
lower classes never exceeds. While below a certain value of k. one needs more number of discs to 
generate required number of lines in higher classes' . The extra discs which are required to fulfill the 
number of lines in the last few classes contribute more for the lower classes (i.e. generaie large number 
of small lines). As one goes to lower classes, the contribution in a lower class from other higher classes 
is more than the act ual strength of that particular class. In order to compensate these excess number of 
lines, a negative value to Ny for that particular class gets generated. 


For an elliptical disc, the probability to find a larger line is lower than that ot a circular disc. 




















If it is possible to assume a certain distribution of disc sizes then the shape and size of the discs 
can be completely determined. For example, figure 6.8 shows that for normally distributed discs, the 
peak in the o/u versus k curves coincide with the actual value of k. Such an assumption may be 
reasonable in view of the fact that many natural processes produce normal distributions. 



rn 


CHAPTER 7 


CONCLUSION 


The major conclusions of this study are summarized below. 

1. Relationships for mean tangent diameter of circular and elliptical discs have been derived. The 
validity of these equations were also tested using simulated data. 

2. The expected length distribution of line segments in the 2-D microstructure has been determined 
from a known size distribution and shape of interfaces (modelled as circular and elliptical discs) 
between particles. 

Conversely, a methodology for transforming the "measured"' line length distribution from observed 
microstructures to size distribution of discs, assuming circular or elliptical shapes. 

4. Without any prior knowledge of shape (circular or elliptical in this case), it has been shown that it 
is possible to extract some limited information on shape of interfaces from the length distribution 
of line segments in the 2-D microstructure. 
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APPENDIX 



APPENDIX A-l 


1. Equation of circle in 3-dimension 

Consider a frame XYZ, which contains a circle of radius R. The center of the circle is at the 
origin of the coordinates. The circle lies on the XY plane and its normal is along Z-axis. Equation of 
circle in the above frame is given by 


X 2 + y 2 = R 2 

and 


( 1 ) 


z = 0 


( 2 ) 


Rotation of coordinates 

Frame Xi YjZj is obtained by rotating frame XYZ by an angle 0 ( clockwise ) around Y-axis. 

Frame X 2 Y 2 Z 2 is obtained by rotating frame XiYiZ, by an angle <j> ( clockwise ) around Z r 
axis. 

So the transformation of coordinates can be written as : 

x = x x cosQ - z x sinQ 

y = y< ( 3 ) 

z = x, sin 0 + z x cosQ 

and 


Xj = x 2 — y 2 sin<b 

= x 2 sin<S> + y 2 cos<b (4) 

7=7 
^ 1 -^2 

From equations (3) and (4), the overall transformation can be written as : 



( 5 ) 


x = c<950 - c<950 sin§ y 2 - sinQ 

y = Sinty X 2 +_a?5(|) y 2 
Z = sinQ cos § xi - sinQ sino y 2 + cos - © 


Substituting the above value of x, y and 2 in equations (1 ) and (2) and writting x->, y 2 , z 2 as x, y, z 
respectively. 


cos 2 Q cos 2 §x 2 - 2cos 2 Q sinfycostyxy - 2 sinQ cosQ cosqxz 
+ cos 2 Q sin® y 2 + wsinQ cosQ sinfyyz + sin 0 ~ 2 
+ sintyx + 2 sin§ costyxy + C os : ® / 

= R 2 


or. 


x 3 (cos 2 Q cos 2 § + sin § ) + y 2 ( cos 2 Q siri® + cos 2 § ) 
+ z sin Q + 2xysin$cos$(l-cos 2 ®) 

+ 2 sinQ cosQ sinfyzy - 2 sinQ cosQ costte 

= R 2 


(7) 


and 


sinQ cos<b x - sinQsinby + cosQ z = 0 


( 8 ) 


If the centre of the circle is at (xo, yo, zo), the equation of circle in 3-dimension can be written as 

(x - xo f. ( cos 2 9 cos 2 <j> + sin <j > ) + (}’ ~ >' 0 / ( cos ' 6 sin 0 + cos 2 * ) 

(z- Za f sin Q + 2(x - x 0 )(y - y 0 )sin§ cosb(l - cos : Q ) (9) 

2(y - y 0 )(z - zo)sinQ cosQ sin<& - 2(x - xo)(z - zo)- sinQ cosQ cosQ 

= R 2 


and 



(x-xo )sinQcos§ - (y - y 0 )sin® sin§ + ( z - Zo )cosQ = 0 


( 10 ) 


2. Expression for end points of the intersecting line 


From equations ( 9) and (10 ), coefficient of various terms can be denoted as : 

cii = cos 2 ® cos 2 § + sin 2 § 
a 2 — cos 2 ® sin § + cos 2 § 
a 3 - (z- zo) sin 9 
aj = sin 2(j) sin © 
as - (z - zo)sin2Q sink 
a 6 - -(z- zo)sin2Q cos§ 

and 

b x = tan 6 

b, = —( z — z 0 )cosQ / sinQ cosfo 

The equations (9) and (10) can be rewritten in the following form : 


( 11 ) 


x' = b,y' + b: 


( 12 ) 


where x' = x-xo ; y' = y-yo and z' = z-zo 

Equation of a plane at height h from the origin and parallel to x-y plane is given by 

Z = h ( 13 ) 

From equations (9), ( 10) and (13), one can write a equation in terms of y as: 


/ '(cubl + ai + aibi) + y (2 cLibibz^ cubi^ cts^ a.6bi) 
= R 2 - as- a,b\~ ctsbi 


(14) 


or 



(15) 


Ay' 2 +By'-C= 0 

where A = aibi 2 + &z + a-ibi , B = 2aibib2 + a4b2 + as +asbi and C = R 2 - as - aib 2 2 - . 


If y{ and yz are the roots of equation (15), the values of roots can be written as : 



-B + -J ~B~ + 4A C 
2A~ 


(16) 


and 


-B-<sIb 2 + 4AC 

2 A 


(17) 


Since ( xo , yo) is the centre of the disc , the y-coordinates of end-points of the intersecting 
line can be written as : 


y, = y'r + y 0 
y 2 = y’r + 


(18) 


Similarity the x-coordinates of end-points of the intersecting line can be written as : 


x, = b, y'r + bi + Xo 

x 2 = b 2 y'v + bz + xo 


(19) 



APPENDIX A-2 


1. Equation of ellipse in 3-dimension 


Consider a Same xyz which contains an elpse whose semi-major and semi-minor axes are 
given by ‘a' and V respectively. The ellipse is oriented in such a manner that its major and mmor axes 
Be along x and y-axes of xyz coordinates respectively. The normal of the elpse lies along 2 -axis and 

the centre of the ellipse is cosidered at the centre of the coordinates. Equation of the ellipse in this 
frame is given by 


and 



( 1 ) 


z - 0 


( 2 ) 


Rotation of Coordinates 

Frame x 1 y J 2 1 is obtained by rotating frame xyz by an angle g (clock wise) around z-axis. 

Frame x^zo is obtained by rotating frame XiyiZi by an angle 0 (clock wise) around y,-axis. 
Frame x 3 y 3 z 3 is obtained by rotating frame X2y 2 z 2 by an angle b (clock wise) around z 2 -axis. 
Transformation of co-ordinates 

x = cos% xi - sinq y t 

y - sin^xi + cos^y, ( 3 ) 

2 = Zi 


and 



xi - cosQ x 2 - sinQ 7 , 

y, = y 2 ( 4 ) 

Zi ~ sinQ X 2 + CO50 7 2 

and 

x 2 = cos§x3-sin§y 3 
y 2 - sin§ x 3 + cos§y 3 (5) 

Z 2 ~ Zs 


From equations (3), (4) and (5), the overall transformation can be written as : 


x - ( cosQ cos (j) cost; - sind sin^Jxs 

-( sint; cos (J) + sinty cosQ cost; )y 3 - sinQ cost; 

y = (sine; cosQ cosfy + cosb ) sin<b )x 3 
+ ( cos<b cosh,- sin^sin^ cosQ )y 3 - sine. sinQ 7 3 


z = sinQ cos® X3 - sinQ sin® y 3 + cosQ (B) 

The equation of ellipse in 3-dimension is obtained by substituting the above values of x. y and z 
in equations ( 1) and ( 2) 



, , f cos 9 coy <j)' cos E, - sin <5> sin c, / ( sin E, coy 6 cos $ + cos E, sin <S) / 

a-:./ ; + 77 i — J 

a b 

, , ( .vm E, co.v (j) + .wztj) co.vG coy E, / f coyE, coyc{) - sinc^si/ub coy 9 f 

+ yj + } 

a b 

2 r cos 2 Z, , sin ^ , . 2Q 

+ zii j — _ + — —} sm y 

a b 

( sin (|) yz'/zE, - cos 9 coy<|) coyE, Jf Ytzz^ coy(|) + Yt'nti) cos 9 cos c , ) 

+ 2 A.» Of, / ; : — 


+ ( yz/z ^ coy 9 cos (j> + cos E, sin (b )( cos E, cos§ - sin E, yz'?z tp coy 9 ) 

b 2 

YZ/Z C, COY cf) + YZ/Z (j) COY 9 COY E ) COS E, YZ/Z E,( YZ7Z E, YZ7Z t|} COY 9 - COY E COY 6 j 

+ 2 y,z.i{ ; + : — : —jsinQ 

a' 1/ 

( Yz'/ztb sint, - coy 9 coy® coyE Jcoyc, ( sin E, coy 9 coy 6 + cos E sin © j yz'tz E, ; . 

+ 2 A - .; z., / ; 7 } «w 

a' b 


yz/z 9 co y(|) xj - 5zh 9 Ym<p ^ + coy9 zj = 0 


2 . Expression for end points of the intersecting line 
Equation (9) can be rewritten as : 

aiX) + Cl2 )'] + O.rzf) + + tZj.Viij + 06 Xj Zj = 1 

where xi = x - xo; yi = y - y«; 7.1 = z - zn; and <11, <12 an arc tlic coefficients of xj , yj .... 

and x?zj respectively. 

Similarily equation (10) can be rewritten as : 


Xj = biy 3 + b 2 


( 12 ) 



where bi = tan <j> and fr> = - zs/cos c}> tan 0 


Equation of a plane in xyz frame which is perpendicular to z-axis and is at a height h from 
the origin, is given by 


Z-h 


From equations (9), (10) and (13 ), one can write a equation in terms of ys as: 


(13) 


(a,b] + <32 + a 4 b,)y 2 , + (2a,bib 2 + cub:* a 5 + a 6 b,) y s = 1 - a 3 - a,bl- a 6 b 2 


or 

Ay) + B y 3 - C = 0 


(14) 


where 


A = a]bj + a 2 + a 4 b] 

B = 2a}blb2 + a4b2 + ay + a^b] 

C = / - a 3 - ajb^- a^b 2 

If yh and y ’2 are the two roots of equation (14) ,the value of roots can be written as: 


y. 


-ff W# : ±4AC 
2 A 


(15) 


v 2 


-B + jB^+JAC 


(16) 


If (xo.yo) is the centre of the disc, the y-coordinates of the end points of the intersecting line 
can be written as : 

>\ = /,+y 0 (17) 

and 

y 2 = y 2 +y 0 (18) 



Similarity the x-coordinates of end points of intersecting line can be written as : 

*1 = + x o 


( 19 ) 


APPENDIX B 


Expressions for di and d 2 defined in section 5.1.1 

General equation of a plane in XYZ coordinates is given by 

Ax + By + Cz + D — 0 (1) 

where the direction of normal to the plane is [A, B, C] and the magnitude of the no rmal is 
given by 

(a 1 + b 2 + or ( 2 ) 

Equation of a plane passing through any point q(Xo, y<>, Zc), whose normal is along the 
direction [A, B, C] is given by 

A(x-xo) + B(y-y 0 ) + C(z- Zo ) = 0 (3) 

or 

Ax + By + Cz - ( Axo + By 0 + Czo ) = 0 (4) 

From equations (1) and (4), one can write equation for D as: 

D = - ( Axo + By 0 + Czo) ( 5 ) 

Figure BJ shows a plane whose normal is given by [A, B, C]. The normal makes an 
angle 9 with Z -axis and the projection of normal on XY plane makes angle <j> with X-axis. 
Consider a unit vector n parallel to the normal [A, B, C]. The equation for the unit vector can 
be written as : 

n — i sin 0 cos (j> + j sin B sin (f> h- k cos 6 
where i , j, k are unit vectors parallel to x, y, z axes respectively. 


( 6 ) 




m 


The vector oq in figure B 1 is given by 

oq = i xo +jy 0 + kzo ( 7 ) 

Magnitude of vector normal to the plane can be written as: 

oq.fi - XoSinOcostp + y 0 sin9sin<f) + Zo cos6 

Let a, p, y be the angles which the normal makes with X, Y and Z axes respectively. 
The direction cosines are given by 

cos a — n.i = sinOcos(f> 
cos ft = n.j = sin 6 sin (p ( 9 ) 

cosy = n.k = cos 6 

From equations ( 1 ) and ( 9 ), the values of A, B, C and D can be written as : 

A = (a 2 + B 2 + C 2 f 2 cos a 

B = (A 2 + B 2 + C 2 ) I/2 cosj3 (10) 

C = (A 2 + B 2 + C 2 f 2 cosy 

The length of normal from the origin to the plane containing first disc is given by 

p, = Xl sinQ, cos (j ) 1 + yposOjSinfi, + zicosd, (11) 

where (xi, y t , zi) and ( 0 i, <j) t ) are the centre and orientation of first disc respectively 

The length of normal from the origin to the plane containing second disc is given by 

p 2 - xisin Q 2 cos(f>, + y,cosQ 2 sin(f > , + z; cos Q 2 (12) 

where (x 2 , y 2 , z 2 ) and (02, <j>2) are the centre and orientation of second disc respectively 

The Equations of planes containing first and second discs can be written as : 

Aix + Biy + Ciz + D, =0 ( 13 ) 


and 


Ai x + Bjy + CiZ + Di ~ 0 


( 14 ) 



where 


Ai ~ Pi sin&jcosfij; 

Bi = Pi sin 0, sin <f> } ; 

Ci = pj cos 0 i 

A2 = P 2 s ™ 62 cos <p 2 ; 

B 2 = P 2 sin O2 sin <f> 2 ; 

C2 = P 2 cos 02 

and 

D, = -(Aixi + B,y, + Cizi) 

D 2 = -(A2X2 + B 2 y 2 + C2Z2 ) 

The equation of the line obtained by intersection of the two planes can be given by the 
combined equation of the planes as follows : 

A,X + B,y + CiZ + D, =_0 

and ( 15 ) 

a 2 x + B 2 y + C 2 z + D2 = 0 

Distances from the centre of discs to the intersecting line can now be written as : 

X 2 + Y 2 + Z 2 
and 

d = ~ ~Y(z 2 -Zo) T + P^( Z 2 - Co) - Z( x 2 - Xo) T + [Y( X 2-Xo) Jr X ( Tr ~ y 0 ) ]' jl/2 

2 x 2 + Y 2 +z 2 

where X = BiC 2 - C1B2 , Y = C1A2 - A t C 2 , Z = A1B2 -B^ and xo, y 0 , zb can be obtained by 
using one of the following eq uati ons 
For X * 0 






APPENDIX C-l , 

This Program Gives the Coordinates of End Points of Lines Obtained 

from Intersection of Discs and Plane 


#include <stdio.h> 

#include <math.h> 

#include <stdlib.h> 

#define pi 3.14159 

/* This program generates the co-ordinates of lines which are stored in 
a file (sample) */ 

int i } j,c,k,n,nv,t; 

float A , B , C , ML , MM , xl , x2 , y 1 , y2 , xO 1 , x02 , yO 1 , y02 ; 
float kk,f 1 ,f2,f3 ,f4,f5,f6,gl ,g2 ,l,h,nl ,randno ; 
float sv , svO , lina , linO , sintheta , dx , dy , dz , rr ,pl , p2 ; 
float al , a2 , b 1 , b2 , c 1 , c2 , dl , d2 , xx , yy , zz , xO , yO , zO , m , p , o , den , dist ; 
float x[1000] ,y [1000] ,z[1000] ,ra[1000] ,rb[1000] ,theta[1000] ,fi[1000] , 
zi [1000] , len [1000] ; 


mainO 

{ 

FILE *f pt ; 

fpt= f open("sample" , "w") ; 

printf ("enter the length of the cube\n") ; 

scanf ("'/.f " ,&1) ; 

printf ("No. of discs/volume\n") ; 
scanf ("7,d" ,&nv) ; 

/* Total no of discs in cube */ 
nl = l*l*l*nv; 
n = nl/1; 

printf ("the value of n= */.d \n",n); 

/* ’n’ will be an integer */ 

printf ("Give ratio b/a\n“); 

scanf ("7.f" ,&kk) ; 

h=0; 

t-0 ; 

c=0 ; 

lin0=0; 

sv0=0; 



for(i=l; i<=n; ++i) 


■C 

randno = (float) rand 0/ (float) RAND_MAX ; 
ra[i] =randno/10.0; 
rb[i]=kk*ra[i] ; 
sv= pi*ra[i]*rb[i] ; 
svO=svO+sv; 

> 

f or (i=l ; i<=n; ++i) 

{ 

found: 

x[i] = 1*( (float) rand 0/ (float) RAND _MAX) ; 
y [i] = 1*( (float) rand()/ (float) RANDJ4AX) ; 
z[i] = 1*( (float) rand ()/ (float) RAND _MAX) ; 

sintheta =((float)RAND_MAX-2 : t : (float)rand())/(float)RAND_MAX ; 

theta [i] = as in (sintheta) +pi/2; 

fi[i] = 2*pi* (float) rand()/ (float) RAND.MAX ; 

zi [i] = 2*pi* (f loat)rand() / (f loat)RAND_MAX ; 

switch (i) 

{ 

case 1 : 
break ; 


default : 

for ( j=l ; j<i; ++j) 

{ 

dx = x[i] - x[j] ; 
dy = y[i] - y [j] ; 
dz = z [i] - z [j] ; 
rr = ra[i] + rb[j] ; 

if (abs(dx)>rr II abs(dy)>rr II abs(dz)>rr) 

i=i; 


else 

{ 

/* distance from the center of the disc to the line */ 

pl=x [i] *sin(theta[i] ) *cos (f i [i] ) + y [i] *cos (theta [i] ) *sin(f i [i] ) 
+ z [i] *cos (theta [i] ) ; 


p2 = x [j] *sin (theta [j] ) *cos (f i [j] ) + y[j]*cos(theta[j])*sin(fi[j]) 
+ z [j] *cos (theta [j] ) ; 


al = pl*sin(theta[i] )*cos(fi[i] ) ; 



a 2 = p2*sin(theta[j] )*cos(fi [j] ) ; 
bl = pl*s in (theta [i] )*sin(f i [i] ) ; 
b2 = p2*sin(theta[j] )*sin(f i [j] ) ; 
cl = pl*cos (theta [i] ) ; 
c2 = p2*cos (theta[j] ) ; 
dl = ~(al*x[i]+bl*y[i]+cl*z[i]) ; 
d2 = ~(a2*x[j]+b2*y[j]+c2*z[j] ) ; 


xx = bl*c2 - cl*b2; 
yy = cl*a2-al*c2; 
zz = al*b2-bl*a2; 


if (xx==0.0) 
if (yy==0.0) 
if (zz==0 . 0) 

printf(" not possible\n") ; 
else 
{ 

x0=(bl*d2-b2*dl) /zz; 
y0=(a2*dl-al*d2) /zz; 
z0=0 . 0 ; 

> 

else 

{ 

x0= (cl*d2-c2*dl) /yy ; 
y0=0 . 0 ; 

z0= (a2*dl-al*d2) /yy ; 

> 

else 

{ 

x0=0 . 0 ; 

y0=(cl*d2-c2*dl)/xx; 
z0= (b2*dl-b l*d2) /xx ; 

> 

/* Distance of the line from the centre of the disc */ 

m=pow( (zz*(y [i] -y0)-yy*(z[i]-z0)) ,2) ; 
p=pow((xx*(z[i]-zO)-zz*(x[i]-xO)) ,2) ; 
o=pow( (yy*(x[i] -x0)-xx*(y [i] -y0) ) ,2) ; 

den = pow(xx,2)+pow(yy,2)+pow(zz,2) ; 


dist=sqrt (m+p+o) / den; 



if (dist>ra[i] ) 

i=i; 

else 

{ 

if (dist<rb[i] H 
++t ; 

goto found; 

> 

else{ 

m=pow ( (zz* (y [ j ] -yO) ~yy* (z [ j ] -zO ) ) , 2) 
p=pow( (xx* (z [j] -zO) -zz* (x [j] -xO) ) , 2) 
o=pow( (yy* (x [j] -xO) -xx* (y [j] -yO) ) , 2) 
dist=sqrt (m+p+o) / den ; 
if (distcrb [j] ){ 

++t ; 

goto found; 

> 

else 

i=i; 

> 


»» 

for(h=0; h<l; h=h+l/10) 

{ 

for(k=l; k<=n; ++k) 

{ 

f l=(pow( (cos (theta [k] )*cos(fi [k] ) *cos (zi [k] )-sin(f i [k] )*sin(zi [k] ) ) ,2) 

/ (ra [k] *ra [k] ) ) + (pow ( (sin (zi [k] ) *cos (theta [k] ) *cos (f i [k] ) +cos (zi [k] ) 
*sin(fi[k])) , 2)/(rb[k]*rb [k])) ; 

f2=(pow( (cos (theta [k] ) *sin(f i [k] )*cos (zi [k] )+cos(f i [k] )*sin(zi [k] ) ) ,2) 

/ (ra [k] *ra [k] ) ) + (pow( (-sin(zi [k] ) *cos (theta[k] ) *sin(f i [k] ) +cos (zi [k] ) 
*cos (f i [k] ) ) , 2) / (rb [k] *rb [k] ) ) ; 

f 3=pow( (h-z [k] ) , 2) *pow (sin (theta [k] ) , 2) *( (pow (cos (zi [k] ) , 2) / (ra[k] *ra[k] ) ) 

+ (pow(sin(zi [k] ) ,2)/(rb [k] *rb [k]))) ; 

f ■ 4=2* ( ( ( (sin(f i [k] ) *sin(zi [k] ) -cos (theta [k] ) *cos (f i [k] ) *cos (zi [k] ) ) * (sin(zi [k] 
*cos (f i [k] ) +sin (f i [k] ) *cos (theta [k] ) *cos (zi [k] ) ) ) / (ra[k] *ra[k] ) ) 

+ ( ( (sin(zi [k] ) *cos (theta [k] ) *cos (f i [k] )+cos (zi [k] ) *sin(f x [k] ) ) *(cos (zi [k] 
*cos (f i [k] ) -sin(zi [k] ) *sin(f i [k] ) *cos (theta [k] ) ) ) / (rb [k] *ro [k] ) ) ) , 

f 5=2*sin (theta [k] ) * (h-z [k] ) * ( ( ( (cos (f i [k] ) *sin(zi [k] ) +cos (thetaCk] ) *sm(f x [k] ) 
*cos (zi [k] ) ) *cos (zi [k] ) ) / (ra [k] *ra [k] ) ) + ( ( (sin (zi [k] ) *cos (thena [k] ) 



*sin (f i [k] ) -cos (zi [k] ) *cos (f i [k] ) ) *sin(zi [k] ) ) / (rb [k] *rb [k] ) ) ) ; 

f 6=2* sin (theta [k] ) * (h-z [k] ) * ( ( ( (sin(f i [k] ) *sin(zi [k] ) -cos (theta [k] ) *cos (f i [k] ) 
*cos (zi [k] ) ) *cos (zi [k] ) ) / (ra [k] *ra[k] ) ) - (( (sin(zi [k] ) *cos (theta [k] ) 

*cos (f i [k] ) +cos (zi [k] ) *sin(f i [k] ) ) *sin(zi [k] ) ) / (rb [k] *rb [k] )) ) ; 

gl=tan(fi [k] ) ; 

g2=- ( (h-z [k] ) *cos (theta [k] ) ) / (sin (theta [k] ) *cos (f i [k] ) ) ; 

A=f l*gl*gl+f 2+f 4*gl ; 

B=2*f l*gl*g2+f ■ 4*g2+f 5+f 6*gl ; 

C=l-f3-fl*g2*g2-f6*g2; 

ML=B*B+4*A*C ; 
if (ML>=0) 

MM=sqrt(ML) ; 
yl= (-B+MM) / (2*A) ; 
y01=yl+y[k] ; 
xl=(gl*yl+g2) ; 
x01=xl+x[k] ; 
y2=-(B+MM)/(2*A) ; 
y02=y2+y [k] ; 
x2=(gl*y2+g2) ; 
x02=x2+x[k] ; 

fprintf (fpt , "*/.f “/.f \n" ,x01 ,y01) ; 

fprintf (fpt , "%f °/.f \n" ,x02 ,y02) ; 

Ien[k]=sqrt(pow((x01-x02) ,2)+pow((y01-y02) ,2)) ; 

IinO=linO+len[k] ; 

++c; 

} 

else 

len [k] =0 ; 

} 

} 

printf ("No . of lines or discs cutting the plane— /(d\n ,c), 

> 



APPENDIX C-2 

This Program Calculates the Length of Various Profile Lines 


#include <stdio.h> 

#include <math.h> 

#include <stdlib.h> 

#define pi 3.14159 

/* This program calculates the length of lines */ 
int i, j , c,k,n,nv,t , q, mm, option, class [50] ; 

float both, lx ,f o ,f ol ,ff , fxl ,fx2, v, mu, sigma, min, max .maxsize ; 
float A,B,C,ML,MM,xl,x2,yl ,y2,x01,x02,y01,y02; 
float kk,f l,f2,f3,f4,f5,f6,gl,g2,l,h,nl,randno; 
float sv , svO , lina, linO , s intheta, dx , dy , dz , rr , pi , p2 ; 
float al , a2 , b 1 , b2 , c 1 , c2 , dl , d2 , xx , yy , zz , xO , yO , zO , m , p , o , den , di st ; 
float x[2000] ,y [2000] ,z[2000] ,ra[2000] ,rb[2000] ,theta[2000] ,fi[2000] , 
zi [2000] ,len [2000] ; 
float gammaCint pp) ; 
main() 

{ 

FILE *f pt ; 

fpt= f openC'out" , "w") ; 

printf ("enter the length of the cube\n") ; 
scanf ("'/.f " ,&1) ; 

printf ("No. of discs/volume\n") ; 
scanf ("7.d" ,&nv) ; 

/* Total no of discs in cube */ 
nl = l*l*l*nv; 
n = nl/1; 

printf ("the value of n= 7.d \n",n); 

/* ’n^ will be an integer */ 

printf ("Give ratio b/a\n") ; 

scanf ("7.f",&kk) ; 

h=0 ; 

t=0 ; 

c=0: 

lin0=0; 

sv0=0 ; _ _ . . . 

printf ("Enter the corresponding # for the type of distnbuxionXn ; ; 



printfC" 1-normal 2-log normal 3-gamma 4-beta 5-weibull \n") 
scant ("’/.d" ,&option) ; 
switch (option) { 
case 1: 

printfC "It takes normal distribution\n") ; 

printf ("the value of mu\n") ; 

scant ("*/.f " ,&mu) ; 

printf ("The value of SIGMA\n") ; 

scant ("*/.f",&sigma) ; 

f or (i=l ; i<=n; ++i) 

{ 

lx=0; 

ff ■ =sigma*sqrt (2*pi) * ( (f loat)rand() / (float) RAND _MAX) ; 
for(fol=0; ff>fol; fol=fol+fo) 

{ 

v=-pow(((lx-mu) /sigma) ,2)/2; 
fxl=exp(v) ; 

lx=lx+((mu+3*sigma)/100) ; 
v=-pow( ( (lx-mu) / sigma) ,2) /2 ; 
fx2=exp(v) ; 

f o= ( (mu+3*sigma) / 100) * (f xl+f x2) *0.5; 

> 

ra[i] =lx; 
rb[i]=kk*ra[i] ; 
sv= pi*ra[i]*rb[i] ; 
sv0=sv0+sv; 

} 

break ; 
case 2: 

printf ("It gives log normal distribution^") ; 
printf ("the value of alf a\n") ; 
scant ("’/.f " ,&mu) ; 
printf ("The value of beta\n"); 
scant ('"/,f" ,&sigma) ; 
for(i=l; i<=n; ++i) 

{ 

lx= (mu+3*sigma) / 100 ; 

f f =sigma*sqrt (2*pi) * ( (float) randO / (float) RAND_HAX) ; 
for(fol=0; ff>fol; fol=fol+fo) 

{ 

v=-pow( ( (log(lx) -mu) /sigma) ,2)/2; 
fxl=exp(v)/lx; 
lx=lx+((mu+3*sigma)/100) ; 
v=-pow(( (log(lx) -mu) /sigma) ,2)/2; 
fx2=exp(v)/lx; 



f 0 = ( (mu+3*sigma) / 100) * (f xl+f x2) *0.5; 

> 

ra[i] =lx; 
rb [i] =kk*ra[i] ; 
sv= pi*ra[i] *rb[i] ; 
sv0=sv0+sv; 

} 

break; 
case 3: 

printfC'It gives gamma distribution^") ; 
printfC'the value of alfa\n") ; 
scanf ("'/.f " ,&mu) ; 
printf ("The value of beta\n") ; 
scanf ("7,f" ,&sigma) ; 
for(i=l; i<=n; ++i) 

{ 

lx=(mu+3*sigma) /100 ; 

ff ■ =pow(sigma ,mu) *gamma(mu) * ( (float) rand() / (f loat)RAND_MAX) ; 
for(fol=0; ff>fol; fol=foi+fo) 

{ 

v=-lx/sigma; 

fxl=exp(v) *pow(lx,mu-l) ; 
lx=lx+((mu+3* sigma)/ 100) ; 
v=-lx/sigma; 

fx2=exp(v) *pow(lx,mu-l) ; 
f o= ( (mu+3*sigma) /100) * (f xl+f x2) *0.5; 

> 

ra[i] =lx; 
rb [i] =kk*ra[i] ; 
sv= pi*ra[i] *rb [i] ; 
sv0=sv0+sv; 

> 

break; 
case 4: 

printfC'It gives beta distribution^") ; 
printfC'the value of alfa\n") ; 
scanf C7.f",&mu) ; 
printf ("The value of beta\n") ; 
scanf ("*/,f" ,&sigma) ; 
for(i=l; i<=n; ++i) 


lx=(mu+3* sigma) /100 ; 

ff =gamma(mu) *gamma(sigma) *( (float) randO / (float )RAND.MAX) /gamma(both) 

f or (f ol=0 ; ff >f ol ; fol=fol+fo) 




v=pow(i-lx J sigma- 1) ; 
fxl=pow(lx,mu-l) ; 
lx=lx+ ( (mu+3*sigma) / 100) ; 

v=pow(l _ lXjSigma-l) ; 
fx2=pow(lx,mu-l) ; 

f o= ( (mu+3*sigma) / 100) * (f xl+f x2) *0.5; 

> 

ra[i] =lx; 

rb[i]=kk*ra[i] ; 

sv= pi*ra[i] *rb[i] ; 
sv0=sv0+sv; 

> 

break; 
case 5: 

printfC'It gives weibull distribution^") ; 
printfC'tbe value of alfa\n“); 
scanf ("‘/.f " ,&mu) ; 
printf ("The value of beta\n") ; 
scanf ("7,f 11 ,&sigma) ; 
for(i=l; i<=n; ++i) 

lx= (mu+3*s igma) / 1000 ; 

ff=((float)rand()/ (f loat)RAND_MAX) / (mu*sigma) ; 
for(fol=0; ff>fol; fol=fol+fo) 

{ 

v=-mu*pow(lx, sigma) ; 
fxl=pow(lx,sigma-l)*exp(v) ; 
lx=lx+( (mu+3*sigma) /100) ; 
v=-mu*pow(lx , sigma) ; 
fx2=pow(lx,sigma-l)*exp(v) ; 
f o=( (mu+3*sigma) /100) *(fxl+f x2)*0 . 5 ; 

> 

ra[i] =lx; 
rb[i]=kk*ra[i] ; 
sv= pi*ra[i] *rb[i] ; 
sv0=sv0+sv; 

> 

break; 
default : 

printf (“Absurd value\n"); 
break;; 

> 

f or (i=l ; i<=n; ++i) 
found: 



X [i] = 1* ((float) rand ()/ (float) RAND _MAX) ; 
y [i] = 1* ( (f loat)rand() / (float) RAND_MAX) ; 
z[i] = 1*( (float) rand ()/ (float) RAND.MAX) ; 

sintheta =((float)RAND_MAX-2*(float)rand())/(float)RAND_MAX ; 
thetaEi] = asin(sintheta) +pi/2 ; 
f i [13 = 2*pi*(f loat)rand()/(f loat)RAND_MAX ; 
zi[i] = 2*pi*(float)rand()/(float)RAND_MAX ; 

switch (i) 

-C 

case 1: 
break; 


default: 

for ( j=l ; j<i; ++j) 

{ 

dx = x[i] - x [ j ] ; 
dy = y[i] - y[j] ; 
dz = z[i] - z[j] ; 
rr = ra[i] + rb[j] ; 

if (abs(dx)>rr II abs(dy)>rr II abs(dz)>rr) 


/* distance from the center of the disc to the line */ 

pl=x [i] *sin(theta[i] )*cos(fi [i] ) + y [i] *cos (theta [i] ) *sin(f i [i] ) 
+ z[i] *cos(theta[i] ) ; 


p2 = x[j]*sin(theta[j])*cos(fi[j]) + y [j]*cos (theta [j])*sin(fi[j]) 
+ z[j]*cos(theta[j]) ; 

al = pl*sin(theta[i] )*cos (f i [i] ) ; 
a2 = p2*sin(theta[j] ) *cos (f i [j] ) ; 
bl = pl*sin(theta[i] )*sin(f i [i] ) ; 
b2 = p2*sin(theta[j] ) *sin(f i [j] ) ; 
cl = pl*cos(thera[i] ) ; 
c2 = p2*cos (thetaCj] ) ; 
dl = ~(al*x[i]+bl*y [i]+cl*z[i] ) ; 
d2 = -(a2*x[j]+b2*y[j] + c2*z[j] ) ; 

xx = bl*c2 - cl*b2; 
yy = cl*a2-al*c2; 
zz = al*b2-bl*a2; 



if (xx==0.0) 
if (yy==0.0) 
if (zz==0.0) 

printfC not possible\n") ; 
else 

x0=(bl*d2-b2*dl)/zz; 
y0= (a2*dl-al*d2) /zz ; 
z0=0 . 0 ; 

> 

else 

{ 

x0=(cl*d2-c2*dl) /yy ; 
y0=0 . 0 ; 

z0= (a2*dl-al*d2) / yy ; 

> 

else 

{ 

x0=0.0; 

y0= (cl*d2-c2*dl) /xx ; 
z0=(b2*dl-bl*d2)/xx; 

> 

/* Distance of the line from the centre of the disc */ 

m=pow ( (zz* (y [i] -yO) -yy* (z [i] -zO) ) ,2) ; 
p=pow((xx*(z[i] -zO)-zz*(x[i]-xO)) ,2) ; 
o=pow( (yy* (x[i] -xO)-xx* (y [i]-yO) ) ,2) ; 

den = pow(xx,2)+pow(yy,2)+pow(zz,2) ; 

dist=sqrt (m+p+o) / den ; 

if (dist>ra[i] ) 
i=i; 
else 
{ 

if (dist<rb [i] ) { 

++t ; 

goto found; 

> * 

else{ 

m=pow((zz*(y[j]-yO)-yy*(z[j]~zO)) ,2) ; 
p=pow((xx*(z[j]-zO)-zz*(x[j]-xO)) ,2) ; 

o=pow( (yy* (x [j] -xO) -xx* (y [j] ~y0) ) ,2) ; 



dist=sqrt (m+p+o)/den; 
if (dist<rb Cj] ){ 

++t ; 

goto found; 

> 

else 

i=i; 

} 

> 

>}}> 

for(h=0; h<l; h=h+l/10) 

{ 

for(k=l; k<=n; ++k) 

{ 

f l=(pow( (cos (theta [k] ) *cos (f i [k] ) *cos (zi [k] ) -sin(f i [k] ) *sin(zi [k] ) ) , 2) / (ra [k] 
*ra[k])) + (pow( (sin(zi [k] ) *cos (theta[k] )*cos (fi [k] )+cos (zi [k] ) 

*sin(fi[k])) ,2)/ (rb [k] *rb [k] ) ) ; 

f2=(pow( (cos (theta [k] ) *sin(f i [k] ) *cos (zi [k] )+cos (f i [k] ) *sin(zi [k] ) ) ,2) / (ra[k] 
*ra [k] ) ) + (pow ( (-sin (zi [k] ) *cos (theta [k] ) *sin (f i [k] ) +cos (zi [k] ) 

*cos (f i [k] ) ) , 2) / (rb [k] *rb [k] ) ) ; 

f3=pow( (h-z Ck] ) , 2) *pow (sin (theta [k] ) , 2) *( (pow (cos (zi [k] ) ,2)/(ra[k] *ra[k] ) ) 

+ (pow(sin(zi [k] ) ,2) / (rb [k] *rb [k] ) ) ) ; 

f 4=2* ( ( ( (sin (f i [k] ) *sin(zi [k] ) -cos (thetaCk] ) *cos (f i [k] ) *cos (zi [k] ) ) * (sin(zi [k] . 
*cos (f i [k] )+sin(f i [k] ) *cos (theta [k] ) *cos (zi [k] ) ) ) / (ra[k] *ra[k] ) ) 

+ ( ( (sin(zi [k] ) *cos (thetaCk] ) *cos (f i [k] )+cos (zi [k] ) *sin(f i [k] ) ) * (cos (zi [k] . 
*cos (f i [k] ) -s in (zi [k] ) *s in (f i [k] ) *cos (theta [kl ) ) ) / (rb [k] *rb [k] ) ) ) ; 

f 5=2*sin (thetaCk] ) * (h-z [k] ) * ( ( ( (cos (f i [k] ) *sin(zi [k] ) +cos (thetaCk] ) *sin(f i [k] ) 
*cos (zi [k] ) ) *cos (zi [k] ) ) / (ra [k] *ra [k] ) ) + ( C <8in(zi W ) (theta [k] ) 

*sin(f i [k] ) -cos (zi [k] ) *cos (f i [k] ) ) *sin(zi [k] ) ) / (rb [k] *rb [k] ) ) ) ; 

f 6=2*sin(theta [k] ) * (h-z [k] ) * ( ( ( (sin (f i [k] ) *sin(zi [k] ) -cos (thetaCk] ) *cos (f : x Ck] ) 

*cos(ziCk]))*cos(ziCk]))/(raCk]*raCk])) - (((sin(zi Ck])*cos (thetaCk]) 

*cos (f i Ck] ) +cos (zi Ck] ) *sin(f i Ck] ) ) *sin(zi Ck] ) ) / (rb Ck] *rb Ck] ) ) ) , 
gl=tan(f i Ck] ) ; 

g2=- ( (h-z Ck] ) *cos (thet a Ck] ) ) / (sin (theta Ck] ) *cos (f i Ck] ) ) ; 

A=f l*gl*gl+f2+f 4*gl ; 

B=2*f I*gl*g2+f4*g2+f5+f6*gl ; 



c=l-f 3-f I*g2*g2-f 6*g2 ; 

ML=B*B+4*A*C; 
if (ML>=0) 

MM=sqrt (ML) ; 
yl=(-B+MM)/(2*A) ; 
yOl=yl+y [k] ; 
x l=(gl*yl+g2) ; 
x01=xl+x[k] ; 
y2=-(B+MM)/(2*A) ; 
y02=y2+y [k] ; 
x2=(gl*y2+g2) ; 
x02=x2+x[k] ; 

fprintf (fpt,'7.f */. f \n" ,x02,y02) ; */ 

len [k] =sqrt (pow ( (x01-x02) , 2) +pow ( (y01-y02) , 2) ) ; 

fprintf (fpt,"7.f \n",len[k]); 

Iin0=lin0+len[k] ; 

++c ; 

> 

else 

len [k] =0; 

> 

> 

printfC'No. of lines or discs cutting the plane= 7,d\n' l ,c); 
Iina=lin0/(1*1) ; 
sv0=sv0/ (1*1*1) ; 
lina=(4*lina) /pi ; 

Iin0=lin0/c ; 

printf (“Average value of line length is= */.f \n" ,lin0) ; 
printf ("4*La/pi = */,f \n“ , lina/10) ; 
printf ("S.area/vol = */.f \n" , svO) ; 

printf ("In how many classes you want to classify the discs\n") ; 
scanf ("7.d" ,&mm) ; 

printf (“Maximum size of discs\n") ; 
scanf ("’/.f" ,&maxsize) ; 
for(q=l; q<=mm; ++q) 
class [q]=0; 
min=0 ; 

max=maxs ize /mm ; 
for(q=l; q<=mm; ++q) 

{ 

f or (i=l ; i<=n; ++i) 

{ 

if (ra[i]>=min && ra[i]<=max) 

++class [q] ; 


else 

i=i; 

> 

min=min+maxs ize/mm ; 
max=max+maxs ize/mm; 

print: f ("No . of discs in '/.d class = */,d\n" ,q, class [q] ) ; 

> 

> 

float gamma(pp) 
int pp; 

int i ; 
float t; 
t=l . 0 ; 
if (pp>l) 

{ 

f or (i=l ; i<=pp-l; ++i) 
t=t*i; 

> 

else 
t=l ; 

return(t) ; 

> 


APPENDIX C-3 

This Program Calculates Number of Intercepts Per Unit Length(NL) on 

a Random Line ° 


#include<stdio . h> 

#include<math . h> 

#include<stdlib ,h> 

#define pi 3.14159 

/* This program is used to find the no of intercepts per unit length(NL) 
This program reads the coordinates from a file (sample) and use them 
to find out the value of ’NL’ */ 
mainO 
{ 

FILE *fpt; 

inh i, j ,k,n,intpt ; 

float c,l, theta, totlen, ml, x0,nl,x[1000] ,y[1000] ,m2[500] ; 

fpt=f openC'sample" , "r") ; 
printf ("Enter the no. of lines\n"); 
scanf ("°/,d" ,&n) ; 
for(i=l; i<=2*n; ++i) 

{ 

f scanf (fpt , "7,f " ,&x [i] ) ; 
fscanf (fpt, "'/.f " ,&y[i]) ; 

> 

printf ("Enter the length of cube\n") ; 

scanf ("*/.f",&l) ; 

totlen=18*l; 

theta=0 ; 

intpt=0 ; 

f or (j = l ; j <=18 ; ++j) 

ml=tan(theta) ; 
for(k=l; k<=n; ++k) 

m2 [k] = (y [2*k] -y [2*k-l] ) / (x [2*k] -x [2*k-l] ) ; 

x0= (m2 [k] *x [2*k-l] -ml*l/2-y [2*k-l] +1/2) / (m2 [k] -ml) ; 

if (x [2*k] >x [2*k-l] ) 



if (x0>=x[2*k-l] && xO<=x[2*k]) 

++intpt ; 
else 

intpt=intpt ; 

> 

else 

if (x0>=x[2*k] && x0<=x[2*k-l]) 

++intpt ; 

else 

intpt=intpt ; 

> 

> 

thet a=t h.eta+( 10/180. 0)*pi; 

> 

nl=intpt/totlen; 

printf ("2*no . of intercept/lengh(nl)= '/.f \n" ,2*nl/10) ; 

/* 2*nl is divided by 10 because in cord2.c we have taken lines 
from 10 sections */ 

} 



APPENDIX C-4 

This Program Obtains the Discrete Distribution of Discs 


#include<stdio . h> 

#include<math . h> 

#include<stdlib .h> 

# define pi 3.14159 

/* This program finds out the distribution of discs from the 
distribution of lines */ 

/* This program calculates for various possible values of b/a 
starting from b/a=l */ 

float dr,kk; 

mainO 

{ 

int i, j ,m,n,naO,tot ,maxno; 
int na[200] ,nv[200] ,np[200] ; 

float nn, sv, coeff ,av, sigma, avg,avgv,sig,sigm,svO,maxlength; 
float min, max, sum, suml, length [80 00] ; 
float kron(int i, int j); 
float pc (int i) ; 

FILE *f pt ; 

FILE *ft ; 

FILE *ft 1 ; 

fpt=f open("out" , "r") ; 

ft=f open("f o" , "a") ; 

ftl=f open("lo" , "w") ; 

printf ("Enter the no. of lines\n"); 

scanf ("7,d" ,&n) ; 

nn=n; 

printf ("The value of n= '/.d\n",n); 
f or (i=l ; i<=n; ++i) 

■C 

f scanf (f pt , "*/,f " , &length [i] ) ; 

} 

printf ("Enter the maxlength of lines\n") ; 
scanf ("7,f " ,&maxlength) ; 



printf ("Enter the no of classes\n"); 
scanf ("7.d" ,&m) ; 
na0=0 ; 
tot=0 ; 

dr=maxlength/m; 
min=0 ; 
max=dr ; 

for(i=l; i<=m; ++i) 

na[i] =0; 

> 

for(i=l; i<=m; ++i) 

f or(j=l ; j<=n; ++j) 

{ 

if (length [j] >=min && length [j] <=max) 

++na[i] ; 
else 
i=i; 

> 

min = min+dr; 
max =max+dr ; 
na[i]=na[i]; 

> 

for(i=l; i<=m; ++i) 

{ 

printf ("na['/.d] =7,d\n" , i,na[i]) ; 

f printf (ft 1 , "7.f \t 7.f\n" ,maxlength*i/10,na[i]/nn) ; 

> 

na0=0 ; 

for(i=l; i<=m; ++i) 
naO =na0+na[i] ; 

f printf (ft , " b/a \t mu \t sigma \t sigma/mu \n\n") ; 
printf ("Total no. of lines=7.d\n" ,naO) ; 

/* naO must be equal to n */ 

/* kron(i,j) is defined as a function*/ 
for(kk=l; kk>=0; kk=kk-0.005) 

{ 

tot=0; 

f or ( j=m; j>0; --j) 
suml=0; 

for(i=m; i>j ; — i) 

{ 

sum = nv[i] *kron(i , j)*pc(i) ; 



suml=suml+sum; 

> 

i=j ; 

nv[i] =(na[j]-suml)/(kron(i, j)*p C (i)) ; 

if (nv[i]>0) 

{ np[i]=nv[i] ; 
nv[i]=nv[i] ; 

> 

else 

{ 

up [i] =nv [i] ; 
nv[i]=0;} 

sv0=0 ; 

for(i=l; i<=m; ++i) 

{ 

sv=nv[i] *pi*kk*i*i*dr*dr/4; 
svO=svO+sv; 

> 

avgv=0 ; 

for(i=l; i<=m; ++i) 

{ 

avg=nv [i] * (2*i-l) *0 . 5*dr*0 . 5 ; 
avgv=avgv+avg ; 

} 

sigm=0 ; 

for(i=l; i<=m; ++i) 

{ 

sig=nv[i]*pow(((2*i-l)*0.5*dr*0.5 - avgv/tot) ,2) ; 
sigm=sigm+sig; 

} 

coeff =sqrt (sigm/tot) / (avgv/tot) ; 
maxno=0 ; 

f or (i=l ; i<=m; ++i) 

{ 

if (np[i] >maxn.o) 
maxno=np [i] ; 
else 

maxno=maxno ; 
s } 

f or (i=l ; i<=m; ++i) 

{ 

if(np[i] < -maxno/10) 

goto last; 

else 

np [i] =np [i] ; 



> 

av=avgv/tot ; 
sigma=sqrt (sigm/tot) ; 

fprintftft,” %1.3f \t\t%f \t tt \t ’/,f \n" ,kk,av, sigma, coeff) 

for(i=l; i<=m; ++i) 

printf ("termimation k=*/.f \n" ,kk) ; 

> 

float kron(i,j) 
int i , j ; 

/* This complete calculates p(i,j)(prob of lines lying in jth 

class after cutting the discs of ith class) */ 

float k,l 3 m,m01 ,m02,th,m0; 

m01=0 ; 

m02=0 ; 

for(th=0; th<=pi/2; th=th+pi/100) 

■C 

k=( (cos (th)*cos (th)+pow( (kk*sin(th) ) , 2) ) * (2* j —1) ) ; 
if ((4*i*i*kk*kk - k*(2*j-l)) >0) 

{ 

l=kk*i*sqrt (4*i*i*kk*kk - k*(2*j-l)); 
m=k/l ; 

> 

else 

m=0 ; 

m01=m01+m; 

> 

for(th=0; th<=pi/2; th=th+pi/2) 

{ 

k=( (cos (th) *cos (th)+pow( (kk*sin(th) ) ,2) )*(2*j-l) ) ; 
if ( (4*i*i*kk*kk - k*(2*j-l)) >0) 

{ 

l=kk*i*sqrt (4*i*i*kk*kk - k*(2*j-l)); 
m=k/l; 

> 

else 

m=0; 

m02=m02+m; 

> 

m0=m01-m02/2; 
m0=m0/50; 
return (mO) ; 

> 

float pc(i) 


{ 

float p,t; 

p=(l-kk*kk)/ (l+kk*kk) ; 

t=4*i*dr*sqrt ( (l+kk*kk) /2) * (l-p*p/ 16-15*p*p*p*p/ 1024) / p 
return (t) ; 

> 



APPENDIX D-l 

Data for Sigma and Mu Values 


The only values of b/a for which any class has more than 107, 
of the discs of the class which has max no of discs are listed 


NORMAL DISTRIBUTION 
mu=l sigma=0.3 

actual b/a=l 
b/a mu 


1.000 

1.031283 

0.995 

1.028964 

0.990 

1 . 026608 

0.985 

1.024167 

0.980 

1.020961 

0.975 

1.017840 

0.970 

1.014654 

0.965 

1.010353 

0.960 

1.005964 

0.955 

1.001018 

0.950 

0.993738 

0.945 

0.983989 

0.940 

0.975867 

0.935 

0.969017 

0.930 

0.960623 

0.925 

0.967177 

0.920 

0.971258 

0.915 

0.963067 

0.910 

0.976605 

0.905 

1.002877 

0.900 

0.998370 

0.895 

1.028552 

0.890 

1.034022 

0.385 

1.028191 

0.880 

1.053045 

0.875 

1.070309 

0.870 

1.074917 


sigma sigma/mu 


0.261989 

0.254042 

0.262205 

0.254824 

0.262403 

0.255602 

0.262655 

0.256457 

0.262854 

0.257457 

0.262960 

0.258351 

0.263029 

0.259231 

0.262797 

0.260104 

0.262817 

0.261259 

0.262255 

0.261988 

0.261284 

0.262931 

0.258871 

0 . 263083 

0.262051 

0.268532 

0.267029 

0.275567 

0 . 273348 

0 . 284553 

0 . 285523 

0.295213 

0 . 297394 

0.306194 

0 . 307800 

0.319604 

0.305130 

0.312440 

0.314780 

0.313877 

0.347921 

0.348489 

0.313246 

0 . 304550 

0.311973 

0.301709 

0.312443 

0.303876 

0.298192 

0.283171 

0.307699 

0 . 287486 

0.293238 

0.272801 



0.865 

1 . 049284 

0.293704 

0 . 279909 

0.860 

1 . 073475 

0.290772 

0 . 270870 

0.855 

1.084835 

0.282114 

0.260053 

0.850 

1.085137 

0.277140 

0.255396 

0.845 

1.085476 

0.272699 

0.251226 

0.840 

1 . 084487 

0.259468 

0.239254 

0.835 

1.071921 

0.273909 

0.255531 

0.830 

1 . 087645 

0.262454 

0.241305 

0.825 

1.100000 

0.262686 

0.238805 

0.820 

1.094122 

0.274197 

0.250610 

0.815 

1.042289 

0.261252 

0.250652 

0.810 

1.093786 

0.273181 

0.249758 

0.805 

1.099341 

0.276005 

0.251064 

0.800 

1 . 104042 

0.278119 

0.251910 

0.795 

1.099872 

0.261770 

0.238001 

0.790 

1 . 103040 

0.285756 

0.259062 

actual 

b/a=.9 



b/a 

mu 

sigma 

sigma/mu 

1.000 

0.968871 

0.252741 

0.260862 

0.995 

0.966860 

0.252941 

0.261611 

0.990 

0.963991 

0.253683 

0.263159 

0.985 

0.961202 

0.253947 

0.264198 

0.980 

0.957972 

0.254713 

0.265887 

0.975 

0.954412 

0.255599 

0.267808 

0.970 

0.950903 

0.255909 

0.269122 

0.965 

0.946410 

0.256502 

0.271026 

0.960 

0.941646 

0 . 257257 

0.273199 

0.955 

0 . 935762 

0.256795 

0.274423 

0.950 

0.929178 

0.256381 

0 . 275922 

0.945 

0.920016 

0.253658 

0.275710 

0.940 

0.912391 

0.256601 

0.281240 

0.935 

0.906130 

0.261366 

0.288443 

0.930 

0.898544 

0.266416 

0.296497 

0.925 

0 . 904028 

0 . 277021 

0 . 306429 

0.920 

0 . 907298 

0.288308 

0.317766 

0.915 

0.899303 

0.297700 

0.331035 

0.910 

0.912196 

0.295950 

0.324437 

0.905 

0.936741 

0.305094 

0.325698 

0.900 

0.931569 

0.335110 

0.359727 

0.895 

0. 960893 

0.305109 

0.317526 

0.890 

0.965760 

0.304252 

0.315039 

0.885 

0 . 959939 

0 . 304389 

0.317092 

0.880 

0.984035 

0.292996 

0.297749 



0.875 

1.000361 

0.303678 

0.303568 

0.870 

1.004831 

0.289547 

0.288154 

0.865 

0 . 980745 

0.288517 

0.294182 

0.860 

1.003552 

0.286486 

0.285472 

0.855 

1.014767 

0.278360 

0.274309 

0.850 

1.015880 

0 . 274297 

0.270009 

0.845 

1.015605 

0.269020 

0.264887 

0.840 

1.015442 

0.256865 

0.252959 

0.835 

1.001990 

0.269749 

0.269213 

0.830 

1.018466 

0.259007 

0.254311 

0.825 

1.031512 

0.258294 

0.250403 

0.820 

1 . 024892 

0.270369 

0 . 263803 

0.815 

0 . 975703 

0.256643 

0 . 263034 

0.810 

1.026060 

0.266720 

0.259946 

0.805 

1.028568 

0.273956 

0.266347 

0.800 

1 . 034924 

0.271550 

0.262386 

0.795 

1.033118 

0.260549 

0.252196 

0.790 

1.035256 

0 . 277759 

0.268299 

0.785 

1.015261 

0.300973 

0.296449 

0.780 

1.045363 

0.273916 

0.262029 

0.775 

1.037518 

0.263029 

0.253518 

actual 

b/a= . 8 



b/ a 

mu 

sigma 

sigma/mu 

1.000 

0.879920 

0.252464 

0.286917 

0.995 

0.878650 

0.252696 

0.287596 

0.990 

0.876706 

0.252605 

0.288130 

0.985 

0.874868 

0.253280 

0.289507 

0.980 

0.872737 

0.253230 

0.290157 

0.975 

0.870708 

0.253162 

0.290754 

0.970 

0.868378 

0.253603 

0.292042 

0.965 

0.866140 

0.253566 

0.292754 

0.960 

0.863455 

0.253624 

0.293732 

0.955 

0.860710 

0.253384 

0 . 294389 

0.950 

0.857633 

0.253227 

0.295263 

0.945 

0 . 853778 

0.252972 

0.296298 

0.940 

0.850257 

0.252836 

0.297364 

0.935 

0.845179 

0.251618 

0.297710 

0.930 

0 . 839642 

0 . 249038 

0.296601 

0.925 

0.835464 

0 . 250441 

0.299762 

0.920 

0 . 828983 

0.250140 

0.301743 

0.915 

0.821286 

0.248136 

0.302131 

0.910 

0 . 822083 

0.255815 

0.311180 

0.905 

0.819707 

0 . 263974 

0.322034 



0.900 0. 

0.895 0. 

0.890 0 . 

0.885 0. 

0.880 0. 
0.875 0, 

0.870 0 

0.865 0 

0.860 0 
0.855 0 

0.850 0 

0.845 0 

0.840 0 

0.835 0 

0.830 C 
0.825 C 
0.820 C 
0.815 C 
0.810 ( 
0.805 < 

0.800 ( 
0.795 1 

0.790 
0.785 1 
0.780 
0.775 
0.770 
0.765 
0.760 
0.755 

actual 

b/a 

1.000 

0.995 

0.990 

0.985 

0.980 

0.975 

0.970 

0.965 

0.960 

0.955 

0.950 


780247 
824518 
828145 
.822386 
.837177 
.811232 
.846788 
. 854457 
.868918 
.887829 
.890256 
.909906 
.921826 
..906392 
1.929658 
(.942608 
) . 922425 
). 925866 
). 939952 
3.927987 
D. 944677 
0.958576 
0.946556 
0.931514 
0.952510 
0.964517 
0.959962 
0.970231 
0.956265 
0 . 976787 


b/a*. 7 
mu 


0.795443 

0.794187 

0.792706 

0.791374 

0.789391 

0.787761 

0.786191 

0.784291 

0.782277 

0.779605 

0.777099 


0.257581 

0.278272 

0.282920 

0.285390 

0.302439 

0.336772 

0.317680 

0.295889 

0.301493 

0.297513 

0.307211 

0.288198 

0.290496 

0.281936 

0.282172 

0.278182 

0.287725 

0.255888 

0.265530 

0.286184 

0.266509 

0.255124 

0.261484 

0.237734 

0.259882 

0.250022 

0.255992 

0.246798 

0.245144 

0.234511 


sigma 

0.256470 

0.256379 

0.256508 

0.256460 

0.256008 

0.256224 

0.256097 

0.256051 

0.255948 

0.255103 

0.254763 


0.330127 

0.337497 

0.341631 

0.347027 

0.361260 

0.415137 

0.375159 

0.346289 

0.346975 

0.335102 

0.345081 

0 . 316734 

0.315131 

0.311053 

0.303522 

0.295120 

0.311922 

0.276377 

0.282493 

0.308392 

0.282117 

0.266149 

0.276248 

0.255213 

0 . 272839 

0.259220 

0.266669 

0.254370 

0.256355 

0.240084 


sigma/ mu 

0.322423 

0.322819 

0.323585 

0.324069 

0.324311 

0.325257 

0.325744 

0.326474 

0 . 327183 

0.327221 

0.327838 



0.945 

0.774410 

0.254666 

0.328852 

0.940 

0.771206 

0.253430 

0.328615 

0.935 

0.767935 

0.252933 

0.329368 

0.930 

0.763581 

0.251222 

0.329005 

0.925 

0.759400 

0.250209 

0.329483 

0.920 

0.753913 

0.248417 

0.329504 

0.915 

0.747254 

0.245206 

0.328143 

0.910 

0.747222 

0.251115 

0.336065 

0.905 

0.743468 

0.254841 

0.342773 

0.900 

0.710563 

0.238761 

0.336017 

0.895 

0.744380 

0.265088 

0.356118 

0.890 

0.746858 

0.270361 

0.361998 

0.885 

0.741239 

0.270697 

0.365195 

0.880 

0.751261 

0.286881 

0.381866 

0.875 

0.721995 

0.307552 

0.425975 

0.870 

0.755999 

0.301357 

0.398621 

0.865 

0.766858 

0.287561 

0.374986 

0.860 

0.778712 

0.294515 

0.378208 

0.855 

0.796348 

0.295318 

0.370840 

0.850 

0.796995 

0.305385 

0.383170 

0.845 

0.817260 

0.294220 

0.360008 j 

0.840 

0.828894 

0.300791 

0.362883 

0.835 

0.813429 

0.292188 

0.359205 

0.830 

0.838558 

0 . 294430 

0.351115 

0.825 

0.851581 

0.289811 

0.340322 

0.820 

0.832200 

0.292414 

0.351375 

0.815 

0.841534 

0 . 274870 

0.326629 

0.810 

0.853630 

0.274923 

0.322063 

0.805 

0.838072 

0.289635 

0.345596 

0.800 

0.859541 

0.275396 

0.320399 

0.795 

0.873851 

0.268160 

0.306871 

0.790 

0.864321 

0.272058 

0.314765 

0.785 

0 . 859267 

0.250010 

0.290958 

0.780 

0 . 868247 

0.263959 

0.304013 

0.775 

0.883290 

0.260408 

0.294816 

0.770 

0 . 878484 

0.263133 

0.299530 | 

0.765 

0.887154 

0.249350 

0.281068 j 

0.760 

0.875896 

0.245708 

0.280522 f 

0.755 

0.897114 

0.237736 

0.265001 i 

.. ' i 

LOG NORMAL 


... ' ■■■ : ' ■ i 

. tj " .... , - ■ A 

■■■' ■ • v 1 


avg=l 

.45 sigma=l . 5 


actual 

b/a=l 



b/a 

mu 

sigma 

sigma/mu 

■ : . . \ ' ' J 

■■■- ' ’ , 1 
■■■ : ■ . f 

. ■ . 1 

i 

■ .-i 



1.000 

1.497808 

1.518613 

1.013890 

0.995 

1.489546 

1.503778 

1.009555 

0.990 

1.487017 

1.504578 

1.011810 

0.985 

1.484466 

1.501763 

1.011652 

0.980 

1.485794 

1.504152 

1.012355 

0.975 

1.483217 

1.504956 

1.014657 

0.970 

1.467933 

1.441786 

0.982188 

0.965 

1.463277 

1.438575 

0.983119 

0.960 

1.453693 

1.420692 

0.977299 

0.955 

1.453846 

1.422221 

0.978247 

0.950 

1.444126 

1.395788 

0.966528 

0.945 

1.436311 

1.390926 

0.968402 

0.940 

1.437482 

1.393150 

0.969160 

0.935 

1.439738 

1.424441 

0.989375 

0.930 

1 . 432749 

1.418933 

0.990358 

0.925 

1.428824 

1.416081 

0.991082 

0.920 

1.417630 

1.385031 

0.977005 

0.915 

1.415648 

1.404859 

0 . 992379 

0.910 

1.416766 

1.421737 

1.003508 

0.905 

1.426316 

1.460801 

1.024178 

0.900 

1.416970 

1.459104 

1.029736 

0.895 

1.428767 

1.491787 

1.044108 1 

0.890 

1.428943 

1.499353 

1.049274 

0.885 

1.415123 

1.477700 

1.044220 

0.880 

1.424961 

1.536637 

1.078371 

0.875 

1.382390 

1.465495 

1.060117 

0.870 

1.422013 

1.559174 

1 . 096456 

0.865 

1.408861 

1.494267 

1.060621 

0.860 

1.426752 

1.569002 

1 . 099702 

0.855 

1.445920 

1.620692 

1 . 120873 

0.850 

1.426010 

1.590120 

1.115084 | 

0.845 

1.456634 

1.647179 

1.130812 

0.840 

1.441176 

1.621003 

1 . 124778 

0.835 

1.417355 

1.598940 

1.128115 

0.830 

1.456694 

1.659685 

1 . 139350 

0.825 

1.482421 

1.684193 

1.136110 

0.820 

1.464489 

1 . 666563 

1.137983 

0.815 

1.444932 

1.606367 

1.111724 

0.810 

1.508784 

1.743770 

1 . 155745 5 

" 1 

0.805 

1.469520 

1.695511 

1.153785 | 

0.800 

1.502058 

1.720663 

1 . 145537 | 

0.795 

1.510017 

1.713334 

1 . 134645 ( 

0.790 

1.507317 

1.731613 

1.148805 | 

0.785 

1.531173 

1.776814 

1 . 160426 

1 

0.780 

1.541975 

1.794839 

1 . 163987 



0.775 

1.523971 

1.763899 

1.157436 

0.770 

1.541516 

1.812773 

1.175967 

0.765 

1.571507 

1.862780 

1 . 185347 

0.760 

1.529630 

1.839608 

1.202649 

0.755 

1.554104 

1.863350 

1.198986 

actual 

b/a 

b/a= .9 

mu 

sigma 

sigma/mu 

1.000 

1.449930 

1.452698 

1.001909 

0.995 

1.443258 

1.443869 

1.000423 

0.990 

1.438505 

1.436752 

0.998782 

0.985 

1.439660 

1.439049 

0.999576 

0.980 

1.434851 

1.435704 

1.000595 

0.975 

1.432000 

1.430306 

0.998817 

0.970 

1.415065 

1.360443 

0.961400 

0.965 

1.411111 

1.357128 

0.961744 

0.960 

1.408116 

1.357221 

0.963856 

0.955 

1.409170 

1.359388 

0.964673 

0.950 

B 

0.945 

1.400000 

1.340398 

0.957427 

. : [ 

1.398972 

1.340398 

0.958130 

0.940 

1.381361 

1.291051 

0.934622 j 

0.935 

1 . 387537 

1.338893 

0.964942 

0.930 

1 . 385373 

1 .339590 

0 . 966953 : 

0.925 

1.383183 

1.353403 

0.978470 ; 

0.920 

1.370393 

1.302728 

0.950624 i 

0.915 

1.365957 

1.328057 

0.972253 i 

0.910 

1.369008 

1.339295 

0 . 978296 

0.905 

1.375269 

1.375645 

1 . 000273 i 

' . ' ! 

0.900 

1.355641 

1.338155 

0.987101 i 

0.895 

1.376087 

1.394350 

1.013271 

0.890 

1.369375 

1.382050 

1.009256 1 

0.885 

1.361417 

1.376956 

1.011414 | 

0.880 

1.380063 

1.458448 

1.056798 

0.875 

1.348562 

1.390637 

1.031200 

0.870 

1.379808 

1 . 484787 

1.076082 | 

0.865 

1.361613 

1.405110 

1.031945 f 

0.860 

1.381818 

1.494425 

1.081492 

0.855 

1.396574 

1.526985 

1.093379 

0.850 

1.390789 

1.531913 

1.101470 

0.345 

1.396529 

1.545456 

1.106641 

0.840 

1.395333 

1.558967 

1.117272 

0.835 

1.365767 

1.507863 

1 . 104041 

0.830 

1.402365 

1.573684 

1 . 122164 



0.825 

1.422504 

1-593435 

1.120162 

0.820 

1.414334 

1-584298 

1.120172 

0.815 

1.397586 

1 - 534722 

1.098124 

0.810 

1.455517 

1-670810 

1.147915 

0.805 

1.425655 

1 • 636619 

1 . 147977 

0.800 

1.439298 

1.618760 

1 • 124687 

0.795 

1.453086 

1-638811 

1.127814 

0.790 

1.457295 

1 - 684435 

1.155864 

0.785 

1.460215 

1.668406 

1.142576 

0.780 

1.475812 

1.712450 

1.160344 

0.775 

1.461538 

1.683421 

1.151814 

0.770 

1.480074 

1.742599 

1 . 177373 

0.765 

1.493680 

1.757754 

1 . 176794 

0.760 

1 . 454459 

1.732394 

1.191092 

0.755 

1.489675 

1.793799 

1.204154 

actual 

b/a=.8 



b/a 

mu 

sigma 

sigma/mu 

1.000 

1.468085 

1.370270 

0.933372 

0.995 

1.467226 

1.371417 

0.934701 

0.990 

1.464586 

1.359138 

0.928002 

0.985 

1.458333 

1 . 348702 

0.924824 

0.980 

1.455633 

1.339525 

0.920236 

0.975 

1 . 452907 

1.337314 

0.920440 

0.970 

1.449261 

1.330651 

0.918158 

0.965 

1 . 446484 

1.328352 

0.918331 

0.960 

1 . 446429 

1.330272 

0.919694 

0.955 

1 . 446372 

1.332207 

0.921068 

0.950 

1 . 443542 

1.329879 

0.921261 

0.945 

1.423965 

1 . 247645 

0.876177 

0.940 

1.423800 

1 . 249327 

0.877460 

0 . 935 

1.418006 

1.237664 

0.872820 

0.930 

1.426292 

1.299968 

0.911432 

0.925 

1.426136 

1.301856 

0.912855 

0.920 

1.409722 

1.257218 . 

0.891819 

0.915 

1.409483 

1.258975 

0.893218 

0.910 

1.412129 

1.277448 

0.904625 

0.905 

1.404070 

1.265841 

0.901551 

0.900 

1.381271 

1.199508 

0.868409 

0.895 

1.392647 

1.250758 

0.898115 

0.890 

1.412099 

1 . 340876 

0.949562 

0.885 

1.403862 

1.330001 

0.947387 

0.880 

1.397611 

1.294688 

0.926358 

0.875 

1.378879 

1 . 305776 

0.946983 



0.870 

1.401817 

1.353235 

0.965344 

0 . 865 

1.392391 

1.288546 

0.925419 

0.860 

1.407343 

1.374898 

0 . 976946 

0.855 

1.411092 

1.386800 

0.982785 

0.850 

1.408688 

1.397255 

0.991884 

0.845 

1.414661 

1.416688 

1.001433 

0.840 

1.421679 

1.440960 

1.013563 

0.835 

1.383167 

1.328498 

0.960475 

0.830 

1.438182 

1.518368 

1.055755 

0.825 

1.438071 

1.487808 

1 . 034586 

0 . 820 

1.436808 

1.518599 

1.056925 

0.815 

1.443182 

1.534049 

1.062963 

0.810 

1.458333 

1.524144 

1.045127 

0.805 

1.451742 

1.555130 

1.071216 

0.800 

1.461648 

1.574972 

1 . 077532 

0.795 

1.471667 

1.565191 

1.063550 

0.790 

1.474038 

1.604396 

1.088435 

0.785 

1.468508 

1.585683 

1.079792 

0.780 

1.486816 

1 . 637742 

1.101509 

0.775 

1.459486 

1.572491 

1.077428 

0.770 

1.495518 

1.712800 

1.145289 

0.765 

1 . 507530 

1.726900 

1.145516 

0.760 

1.451176 

1.625682 

1.120252 

0.755 

1.463144 

1.625058 

1.110661 

0.750 

1.365761 

1.736296 

1.271303 

0.745 

1.499739 

1.733397 

1.155799 

0.740 

1.471928 

1.735260 

1.178903 

0.735 

1.533996 

1.694605 

1 . 104700 

0.730 

1.505901 

1.789793 

1.188519 

0.725 

1.605084 

1.805933 

1.125133 

0.720 

1 . 579449 

1.741669 

1 . 102707 

0.715 

1 . 633457 

1.774163 

1.086140 

0.710 

1.621269 

1 . 745726 

1.076766 

0.705 

1.569568 

1.708683 

1.088633 

actual 

b/a= .7 



b/a 

mu 

sigma 

sigma/mu 

1.000 

1.450079 

1.338959 

0.923370 

0.995 

1.449132 

1.340062 

0.924734 

0.990 

1.446315 

1.337791 

0.924965 

0.985 

1.445350 

1.338897 

0 . 926348 

0.980 

1.445288 

1.340886 

0 . 927764 

0.975 

1.447115 

1.342644 

0 . 927807 

0.970 

1.438607 

1.331664 

0.925662 



0.965 

1.440428 

1.333428 

0.925716 

0.960 

1.427032 

1.295056 

0.907517 

0.955 

1.417414 

1.246707 

0.879564 

0.950 

1.417213 

1.248426 

0.880902 

0.945 

1.418997 

1.250089 

0.880967 

0.940 

1.411010 

1.238125 

0 . 877474 

0.935 

1.404950 

1.221658 

0.869538 

0.930 

1.395903 

1.208155 

0.865501 

0.925 

1.404412 

1.258157 

0.895860 

0.920 

1.404139 

1.259967 

0 . 897324 

0.915 

1.397920 

1.256015 

0.898488 

0.910 

1.391638 

1 . 243602 

0.893624 

0.905 

1.380155 

1.217082 

0.881844 

0.900 

1.365484 

1 . 175053 

0.860539 

0.895 

1.374130 

1.206324 

0.877881 

0.890 

1.370629 

1.215873 

0.887091 

0.885 

1.370167 

1.217584 

0.888639 

0.880 

1.376770 

1 . 252783 

0 . 909944 

0.875 

1.343750 

1.196593 

0.890488 

0.870 

1.373656 

1.299512 

0 . 946024 

0.865 

1.366892 

1.238888 

0.906354 

0.860 

1.382246 

1.331643 

0 . 963390 

0.855 

1.385949 

1.344127 

0 . 969824 

0.850 

1.385550 

1.354741 

0.977764 

0.845 

1.398063 

1.389262 

0.993705 

0.840 

1.398699 

1.401016 

1.001656 

0.835 

1.361842 

1 . 284968 

0.943551 

0.830 

1.405603 

1.421148 

1.011060 

0.825 

1.420841 

1.446723 

1.018216 

0.820 

1.406011 

1.419472 

1.009574 

0.815 

1.419146 

1.457975 

1 . 027361 

0.810 

1.434730 

1.471403 

1 . 025560 

0.805 

1.429961 

1.519097 

1 . 062334 

j 

0.800 

1.436644 

1.483806 

1.032828 

0.795 

1.460630 

1 . 539442 

1.053958 | 

0.790 

1 . 459493 

1 . 575886 

1 . 079749 

0.785 

1.470000 

1.588938 

1.080910 

0.780 

1.461869 

1.550018 

1.060299 

0.775 

1.447597 

1.546021 

1.067991 

0.770 

1.476337 

1.641521 

1.111888 

0.765 

1.484959 

1.659514 

1.117549 i 

' ’ ’ , „• * , "v n 1 

0.760 

1.457105 

1.632675 

1 . 120492 j 

0.755 

1.459574 

1.615550 

1 . 106864 

0.750 

1.353700 

1.666857 

1.231335 

0.745 

1.478448 

1.666599 

1 . 127262 



0.740 

1.460886 

1.677217 

1.148082 

0.735 

1.532667 

1.695497 

1.106239 

0.730 

1.488082 

1.729798 

1 . 162435 

0.725 

1.594697 

1.764826 

1.106684 

0.720 

1.583424 

1.759755 

1.111361 

0.715 

1.610230 

1.708015 

1.060727 

0.710 

1.605132 

1.693349 

1.054959 

0.705 

1.570205 

1 . 694949 

1 . 079444 

GAMMA DISTRIBUTION 


actual 

b/a=l 

AVG=1.6 sigma=l 

.2 

b/a 

mu 

sigma 

sigma/mu 


1.000 

1.712597 

1.307424 

0.763416 

0.995 

1.710124 

1.305037 

0.763124 

0.990 

1.705334 

1.304846 

0.765156 

0.985 

1.698834 

1.298396 

0.764287 

0.980 

1.695096 

1.296614 

0.764921 

0.975 

1.691318 

1.295586 

0.766021 

0.970 

1.682827 

1.285583 

0.763942 

0.965 

1 . 672863 

1.275157 

0.762260 

0.960 

1.669595 

1.274811 

0.763545 

0.955 

1.663690 

1 . 268447 

0.762429 

0.950 

1.655308 

1.263418 

0.763252 

0.945 

1 . 640776 

1.244226 

0.758315 

0.940 

1 . 628472 

1.232318 

0 . 756733 

0.935 

1.618873 

1.222821 

0.755353 

0.930 

1 . 607458 

1.214729 

0.755683 

0.925 

1.604167 

1.220945 

0.761108 

0 . 920 

1.599159 

1.226759 

0.767127 

0.915 

1.579299 

1.207875 

0.764817 

0.910 

1.584271 

1.225495 

0.773539 

0.905 

1 . 597354 

1.262653 

0 . 790465 

0.900 

1 . 565278 

1.254129 

0.801218 

0.895 

1.606324 

1.300264 

0.809466 

0.890 

1.601993 

1 . 306750 

0 . 815703 

0.885 

1.591447 

. 1.305536 

0.820345 

0.880 

1.600189 

1.331216 

0.831911 

0.875 

1.553363 

1.340118 

0.862720 

0.870 

1.609110 

1 . 375339 

0.854720 

0.865 

1.610027 

1.344321 

0.834968 

0.860 

1 . 631298 

1.382157 

0 . 847274 

0.855 

1 . 666570 

1.417771 

0.850712 

0 . 850 

1 . 649038 

1.415390 

0.858312 



0.845 

1.685541 

1.437537 

0.852864 

0.840 

1.687698 

1.445501 

0.856493 

0.835 

1.652621 

1.426269 

0.863034 

0.830 

1.706041 

1.453931 

0.852225 

0.825 

1.747231 

1.467623 

0.839971 

0.820 

1.712018 

1.445000 

0.844033 

0.815 

1.692565 

1.389355 

0.820858 

0.810 

1.768840 

1.451953 

0.820850 

0.805 

1.731413 

1.458644 

0.842458 

0.800 

1.787460 

1.478620 

0.827219 

0.795 

1.815731 

1.483481 

0.817016 

0.790 

1.802900 

1.479101 

0.820401 

0.785 

1.806567 

1.464752 

0.810793 

0.780 

1 . 844672 

1.497783 

0.811951 

0.775 

1 . 839837 

1.506424 

0.818781 

0.770 

1.852768 

1.517568 

0.819081 

0.765 

1.916701 

1.554318 

0.810934 

0.760 

1.874247 

1.561174 

0.832961 

0.755 

1.923801 

1.562799 

0.812350 

0.750 

1 . 904756 

1.709543 

0.897512 

0.745 

1.928017 

1 . 585772 

0 . 822488 

0.740 

1.932162 

1 . 644921 

0.851337 

0.735 

1.971341 

1.569683 

0.796251 

0.730 

1.905242 

1.600071 

0.839826 

0.725 

1 . 979305 

1.543024 

0.779578 

0.720 

1 . 977972 

1.586710 

0.802190 

0.715 

2.005569 

1 . 547634 

0.771668 

0.710 

1.980071 

1.530373 

0.772888 

actual 

CD 

II 



b/a 

mu 

sigma 

sigma/mu 

1.000 

1.665432 

1.203792 

0.722811 

0.995 

1 . 659776 

1.196039 

0.720603 

0.990 

1.656576 

1.195330 

0.721567 

0.985 

1.649215 

1.188217 

0.720475 

0.980 

1 . 648684 

1 . 189202 

0.721303 

0.975 

1.642361 

1 . 185422 

0.721779 

0.970 

1.637134 

1 . 184249 

0.723367 

0.965 

1.629512 

1 . 176730 

0.722136 

0.960 

1 . 624579 

1.172655 

0.721821 

0.955 

1.619164 

1.171297 

0.723396 

0.950 

1.611310 

1.163366 

0.722001 

0.945 

1.599366 

1.145154 

0.716005 



o . 940 

1.594095 

1.144199 

0.717774 

O . 935 

1.579016 

1.128659 

0.714786 

0 .930 

1.573529 

1.129363 

0.717726 

0 .925 

1 . 570434 

1.131399 

0 . 720437 

0 .920 

1.556455 

1.121729 

0.720695 

O .915 

1.540979 

1 . 105445 

0.717365 

O .910 

1.540101 

1.118718 

0.726393 

O .905 

1.550327 

1.151817 

0.742951 

O .900 

1.517699 

1 . 131474 

0.745519 

O .895 

1.551254 

1 . 178793 

0.759897 

0 .890 

1.542969 

1.178188 

0.763585 

0 .885 

1.531907 

1.181902 

0.771523 

O .880 

1.540185 

1.206749 

0.783509 

0 .875 

1.478777 

1.198674 

0.810584 

0.870 

1.535301 

1.235997 

0.805052 

0 .865 

1 . 542407 

1.215892 

0.788308 

0.860 

1.563469 

1.252661 

0.801207 

0.855 

1.593314 

1 . 284767 

0.806349 

0.850 

1.575689 

1 . 284082 

0.814934 

0.845 

1.612853 

1.308070 

0.811029 

0.840 

1.616185 

1.324116 

0.819285 

0.835 

1 . 576838 

1.291485 

0.819035 

0.830 

1.639055 

1.335520 

0.814811 

0.825 

1.676622 

1.346917 

0 . 803352 

0.820 

1 . 644494 

1.325080 

0.805768 

0.815 

1.634421 

1.280550 

0.783489 

0.810 

1.712381 

1 . 340458 

0.782804 

0.805 

1 . 654767 

1.325214 

0.800846 

0.800 

1.726008 

1.362148 

0.789190 

0.795 

1.756341 

1.367685 

0.778713 

0.790 

1.744300 

1.366481 

0 . 783398 J 

0.785 

1.759248 

1.348030 

0.766253 | 

0.780 

1.774306 

1.371019 

0.772708 

0.775 

1.784884 

1.391993 

0.779879 

0.770 

1.783654 

1.389782 

0.779177 

0.765 

1.845046 

1.419792 

0.769516 

0.760 

1.799613 

1.421313 

0.789788 : 

0.755 

1.856015 

1.432171 

0.771638 j 

actual 

b/a 

b/a= . 8 

mu 

sigma 

y ' ' " y f 

sigma/mu 

1.000 

1 . 555249 

1 . 170462 

0 . 752588 

0.995 

1 . 549826 

1.161369 

0.749355 

0.990 

1.542118 

1 . 153376 

0.747917 



°- 985 1.541053 

1.154109 

0.748909 

°- 980 1.537725 

1.153581 

0.750187 

°- 975 1.534076 

1.151461 

0.750589 

°- 970 1.528690 

1 . 146847 

0.750216 

0 • 968 1.520670 

1.138314 

0.748561 

°- 960 1.516283 

1.136338 

0 . 749424 

°- 955 1.509259 

1.131310 

0.749580 

0.95° 1.506815 

1.132582 

0.751639 

0.945 1.495476 

1.117767 

0.747432 

0.940 1.488208 

1.112321 

0.747423 

0.935 1.481089 

1 . 102773 

0.744569 

0.930 1.468109 

1.087360 

0.740654 

0.925 1.461743 

1.085492 

0.742601 

0.920 1.452063 

1.077165 

0.741817 

0.9 1 5 1.439773 

1.067834 

0.741669 

0.9 10 1.445411 

1.090129 

0 . 754200 

0.905 1.456064 

1.123195 

0.771391 

0.900 1.433923 

1 . 109337 

0.773638 

0.89 5 1.453305 

1.142198 

0.785932 

0.890 1.450325 

1.146390 

0.790436 

0.88 5 1.435128 

1 . 138922 

0.793603 

0.880 1.442655 

1.168051 

0 . 809654 

0.875 1.386518 

1.152892 

0.831502 

0.87 0 1.439749 

1.201142 

0 . 834272 

0.855 1.440132 

1.171899 

0.813745 

0.860 1.459214 

1.213078 

0.831323 

0.85 5 1.485911 

1 . 243524 

0 . 836876 

0.85 0 1.470859 

1.246252 

0 . 847295 

0.845 1.496076 

1.259161 

0.841642 

0.840 1.496711 

1.274879 

0.851787 

0.835 1.460173 

1.237896 

0 . 847773 

0.830 1.517578 

1.287118 

0.848140 

0.825 1.558043 

1.310218 

0 . 840938 

0.820 1.524371 

1.283617 

0.842063 

0.815 1.504471 

1.222455 

0.812548 

0.810 1.589606 

1 . 307239 

0.822367 

0.805 1.525983 

1 . 276022 

0.836196 

0.800 1.600931 

1 . 324319 

0.827218 

0.795 1.623118 

1 . 322966 

0.815077 

0.790 1.609776 

1.319946 

0.819956 

0.785 1.633100 

1 . 313827 

0 . 804499 

0.780 1.638508 

1.328368 

0.810718 

0.775 1.669352 

1.376955 

0 . 824844 




actual b / a =.7 



b/ a 

mu 

sigma 

sigma/mu 

1.000 

1.457544 

1.026189 

0.704054 

0.995 

1.457541 

1.026929 

0 . 704563 

0.990 

1.457538 

1.027674 

0.705076 

0.985 

1 . 452740 

1.022956 

0.704157 

0.980 

1.452717 

1.023694 

0.704676 

0.975 

1.451897 

1.023955 

0.705253 

0.970 

1 . 449446 

1.019842 

0.703608 

0.965 

1 . 447772 

1.020343 

0.704768 

0.960 

1 . 447727 

1.021084 

0.705301 

0.955 

1.442767 

1.016139 

0.704299 

0.950 

1.442701 

1.016871 

0 . 704839 

0.945 

1.442635 

1.017609 

0.705382 

0.940 

1.435897 

1.003615 

0.698946 

0.935 

1.431590 

0.998846 

0.697718 

0.930 

1.431475 

0.999529 

0.698251 

0.925 

1.426301 

0.994151 

0.697014 

0.920 

1.424419 

0.994494 

0.698175 

0.915 

1.424270 

0.995171 

0 . 698723 

0.910 

1.419787 

0.990138 

0.697385 

0.905 

1.417836 

0.990451 

0.698565 

0.900 

1.413262 

0.985276 

0.697165 

0.895 

1.406016 

0.969678 

0.689663 

0.890 

1.401316 

0.964142 

0.688027 

0.885 

1.399206 

0.964273 

0.689157 

0.880 

1.395198 

0.959087 

0 . 687420 

0.875 

1 . 380376 

0.924559 

0.669788 

0.870 

1.380796 

0.935655 

0 . 677620 

0.865 

1.381221 

0 . 946755 

0 . 685448 

0 . 860 

1.376174 

0 . 940575 

0.683471 

0.855 

1.373819 

0.945617 

0.688313 

0.850 

1 . 367448 

0.939253 

0.686866 

0.845 

1.367800 

0 . 950657 

0.695026 

0.840 

1.357661 

0 . 937246 

0.690338 

0.835 

1.342427 

0.915927 

0.682292 

0.830 

1.351230 

0.936864 

0.693342 

0.825 

1.361431 

0.957528 

0.703325 

0.820 

1 . 362562 

0.959293 

0.704036 

0.815 

1.357859 

0 . 964037 

0.709968 

0.810 

1.380034 

1.009589 

0.731568 

0.805 

1 . 364407 

0.979555 

0.717935 

0.800 

1 . 376706 

1.016815 

0.738585 

0.795 

1.384036 

1.041638 

0.752609 

0.790 

1.379340 

1.036750 

0.751628 



f 


0.785 

1.385927 2T ICU 

1.060964 

0.765527 

0.780 

1.387346 

1.063494 

0.766568 

0.775 

1.364247 

1.048711 

0.768710 

0.770 

1.365506 

1.056507 

0.773711 

0.765 

1.381831 

1.105160 

0.799780 

0.760 

1.340149 

1.058936 

0.790163 

0.755 

1.348684 

1.102165 

0.817215 

0.750 

1.158401 

0.992965 

0.857186 

0.745 

1.369782 

1.106988 

0.808150 

0.740 

1 .336946 

1.066716 

0.797875 

0.735 

1.439623 

1.147371 

0.796994 

0.730 

1.359406 

1.066959 

0 . 784872 

0.725 

1.511961 

1.182563 

0.782139 

0.720 

1.498106 

1.152647 

0.769403 

0.715 

1.551887 

1.186683 

0.764671 

0.710 

1.566399 

1.199250 

0.765610 

0.705 

1.478206 

1.225510 

0.829052 

0.700 

1.557261 

1.206976 

0.775063 

0.695 

1.587209 

1.209630 

0.762111 

0.690 

1.659483 

1.247126 

0.751515 

0.685 

1.585500 

1 . 238876 

0.781379 

0.680 

1.617110 

1.178241 

0.728609 

0.675 

1.722868 

1.267939 

0.735947 

0.670 

1.671092 

1.257694 

0.752618 

0.665 

1.705040 

1.215641 

0.712969 

0.660 

1.709012 

1.279517 

0.748688 

0.655 

1.770916 

1.272506 

0.718558 

0.650 

1.700884 

1 . 256257 

0.738591 

0.645 

1.817036 

1.284618 

0.706985 

0.640 

1.770094 

1.278570 

0.722318 

0.635 

1.835714 

1.238256 

0.674536 

0.S30 

1.764768 

1.209625 

0.685430 
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